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Abstract 

In this article we continue our study of chiral fermions on a quantum curve. 
This system is embedded in string theory as an I-brane configuration, which 
consists of D4 and D6-branes intersecting along a holomorphic curve in a 
complex surface, together with a 5-field. Mathematically, it is described by a 
holonomic P-module. Here we focus on spectral curves, which play a promi- 
nent role in the theory of (quantum) integrable hierarchies. We show how 
to associate a quantum state to the I-brane system, and subsequently how 
to compute quantum invariants. As a first example, this yields an insightful 
formulation of (double scaled as well as general Hermitian) matrix models. 
Secondly, we formulate c = 1 string theory in this language. Finally, our for- 
malism elegantly reconstructs the complete dual Nekrasov-Okounkov partition 
function from a quantum Seiberg-Witten curve. 
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1 Introduction 

The study of two-dimensional conformal field theories on Riemann surfaces has had 
many fruitful applications in physics and mathematics. In many ways the field theory 
of a chiral free fermion is the most important and instructive example. In that case 
one considers a Riemann surface or algebraic curve E together with a line bundle C 
that comes equipped with connection A. The free fermion partition function computes 
the determinant of the twisted Dirac operator dA coupled to the line bundle C This 
determinant has many interesting properties, e.g. the dependence of this determinant on 
the connection A is captured by the Jacobi thcta-function. It is known for a long time 
that these chiral determinants arc closely related to integrablc hierarchies of KP-type 
[1, 2, 3]. In the simplest case this relation arises as follows. One picks a point P G S on 
the curve together with a local trivialization e* of the line bundle around P. The ratio 
with respect to a reference connection Aq 

det dA 

T{t) — = — 

det 

then becomes a so-called tau-function of the KP-hierarchy. In the Hamiltonian formula- 
tion one associates a state | W) in the fermionic Fock space T to the hne bundle on E — P. 
In the semi-infinite wedge representation of the Fock space this state can be considered 
as the wedge product of a basis that spans the space of holomorphic sections 

W = //°(E-P,£). 

Similarly, a coherent state \t) is associated to the local trivialisation around P. Combining 
these two ingredients the tau-function can be written as 

T{t) = (1.1) 

With the advent of matrix models it became clear that string theory can also give rise 
to solutions of KP-type of the form (1.1). More recently this connection to integrable 
hierarchies has been reformulated and generalized through the methods of topological 
strings [4, 5]. These string theory solutions are similar, but not equivalent, to the familiar 
geometric solutions coming from CFT that are sketched above. In particular the relevant 
Fock space state |>V) does not have a purely geometric interpretation as generated by a 
space of sections over a curve. Yet, in the string theory setting an algebraic curve E does 
appear. (Here it should be stressed that this curve is not a string world-sheet, but should 
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be considered as (part of) the target space geometry.) But in this case there is an extra 
parameter: the string couphng constant A. Only in the genus zero or classical limit A — > 
a geometric curve arises. There have been many indications that A should be interpreted 
as some form of non-commutative deformation of the underlying algebraic curve. In the 
simplest cases E appears as an affine rational curve given by a relation of the form 

F(x,|/)=0, 

in the complex two- plane C^, with a (local) parametrization 

with p, q polynomials. Of course, p and q commute: [p, q] = 0. However, the string-type 
solutions with A 7^ are characterized by quantities P and Q that no longer commute 
but instead satisfy the canonical commutation relation 

[P, Q] = A. 

In this case clearly P, Q cannot be polynomials, but are represented as differential oper- 
ators, i.e. polynomials in z and dz. As we will point out in this paper a suitable concept 
to frame these solutions is a "D-module. Instead of classical curve in the {x, y)-plane, we 
should think of a quantum curve as an analogue in the non-commutative plane [x, y] = A. 
If we interpret 

one can identify such a quantum curve as a holonomic "D-module W for the algebra 
V of differential operators in x. Now there is a straightforward way in which such a 
"D-module gives rise to a solution of the KP-hierarchy. By definition W carries an action 
of both X and d^- However we are free to ignore the second action, which leaves us 
with the structure of an O-module, O being the algebra of functions in x. By applying 
the infinite-wedge construction to the module W we obtain in the usual way a state 
\yV) in the fermion Fock space. Roughly speaking, W can be considered as the space 
of local sections that can be continued as sections of a (non-commutative) P-module, 
instead of sections of a line bundle over a curve. This set-up can be generalized in 
many ways and in this fashion several constructions in topological string theory, matrix 
models and integrable hierarchies can be connected. It is the purpose of this paper to 
explain the connections between these familiar ingredients from the "D-module perspective. 
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This paper is structured as follows: 

In section 2 we introduce our notion of a quantum curve and provide a construction 
of a tau-function associated to it. This tau-function arises, in an appropriate sense, from 
a quantization of the Krichever correspondence described in § 2.1. The physical system 
relevant for this quantization consists of an intersecting brane configuration with S-field 
in string theory. It is introduced in [5] and reviewed in § 2.2. The D4 and D6-branes wrap 
an affine complex curve E that is embedded in a complex symplectic plane. Endpoints of 
the strings stretched between D4 and D6-branes appear as fermionic modes on E, which 
are quantized by the B-Held. This turns the chiral fermions into sections of a so-called 
"D-module. In § 2.2 we explain in which sense a "D-module quantizes the spectral curve and 
in § 2.3 we discuss how one can associate a fermionic state |W), and thus a tau-function, 
to such a "D-module. 

In sections 3, 4 and 5 we analyse three physical systems in which above quantization 
is realized: respectively matrix models, c — 1 string theory and H — 2 supersymmetric 
gauge theories. We find that a quantization of the underlying classical curve yields a 
differential system that determines the corresponding partition functions. In other words, 
we see how our formalism in section 2 gives a unifying picture of these topics in terms of 
a underlying quantum curve. 

2 Quantum curves and invariants 

The main object of interest in this paper is a chiral fermion field living on a holomorphic 
quantum curve. This set-up is embedded in string theory as a configuration of D4 and 
D6-branes that intersect along a classical curve E. Turning on a B- field on the D6- 
brane quantizes the curve E. As was shown in [5] this intersecting brane configuration 
is closely related to topological string theory, supersymmetric gauge theory and matrix 
models. More precisely, it relates to the topological B-model on non-compact Calabi-Yau 
backgrounds of the form 

Xy, : uv — F{z, w) — 0, 

that is modeled on an afiine curve S defined by the equation F{x,y) = 0. The topological 
string partition function admits an expansion 
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in the topological coupling constant A, whose classical contribution captures the com- 
plex periods 

JAi JBi 

of Xy,, while the semi-classical contribution is known to compute a chiral determinant 

exp^i = det9s (2.1) 

on E. All higher order Tg's give quantum corrections to these results. 

As we alluded to in the introduction, the chiral determinant (2.1) has an elegant 
interpretation in terms of certain geometric solutions of the KP hierarchy, which is known 
as the Krichever correspondence. In this context the chiral determinant is known as a 
tau-function. On the other hand, the total topological string partition function is also 
known to represent a tau-function of a KP hierarchy, though in this case it doesn't have a 
similar geometric interpretation. The aim of this section is to propose a quantum analog 
of the Krichever correspondence, starting from a quantum curve. We conjecture that this 
prescription computes the all-genus topological string partition function. 

In this section we start by reviewing the Krichever correspondence. We continue by 
reviewing the intersecting brane system and explain what we mean by a quantum curve. In 
the last subsection we line out our prescription to obtain invariants from such a quantum 
curve. 

2.1 Krichever correspondence 

In this section we review the geometric Krichever correspondence that underlies the 
genus 1 free energy Ti of the topological string. In the simplest scenario we start with a 
Riemann surface E with a single puncture P. We study a chiral fermion field 

reZ+1/2 

on E which is coupled to a line bundle C. The Hilbert space V, of this fermion field is built 
by acting on the Dirac vacuum |0) with the fermionic modes and ■0]!. The boundary 
conditions near P, i. e. the choice of local coordinates and a choice of local frame of C 
near P, are encoded in a coherent state 

{t\ = (0|e^*""", 



6 



where q;„ = : : are the bosonized modes. The partition function of the fermion 

field sweeps out a state | W) in the Hilbert space V,, and for a given choice t of boundary 
conditions it reads 

Tit) = {t\W). (2.2) 

The Krichever correspondence tells us precisely how to find the state |>V). Choosing 
as a local coordinate around P, we define a subspace 

W = H\J: - P, jC) G C[z] ® C[[z-% (2.3) 

By picking a semi-infinite basis Wn = z"'{l + 0{z~^)) of this subspace, it can be quantized 
into a fermionic state 

\W) ^WiAw2Aw3A...en. (2.4) 
on which the fermionic modes act as 

Any state |>V) that we obtain in this way looks hke 

I W) = ^|0), where g = exp Cnm^ni^l) e Gl{(x,). (2.5) 

All states of the above form parametrize an infinite Grassmannian, which is well-known 
to give an elegant geometric formulation of the KP integrable hierarchy. (A more detailed 
review of these issues can be found in Appendix A, which will be useful later on.) 

Important for now is that although one can associate a tau function to any element 
I W) in the Grassmannian, as in equation (2.2), only a dense subset of subspaccs W in the 
infinite Grassmannian of measure allows for a geometric Krichever interpretation. This 
subset can be characterized as follows. Basically, a subspace W has a geometric origin 
when there is an algebra A such that 

A-WCW, (2.6) 

with A being non-trivial, i.e. A ^ C. In this situation the underlying curve E can be 
defined in terms of its spectrum A = i7°(E — P,C) and the tau-function (2.2) has an 
interpretation as a fermionic determinant det d^. 
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2.2 B-field, "D-modules and quantum curve 

Our motivation for quantizing the Krichever correspondence comes from [5] , where it 
is argued that topological string theory is related by a chain of duahties to a system of 
intersecting D4 and D6 branes in type IIA with a background B-&e\d. We now review 
some aspects of this relation. 

I-brane configuration 

In the intersecting brane set-up a central role is played by a holomorphically embedded 
curve E C B, given by the equation 

E : F{z,w) = 0, 

where B = C x C (or possibly with either C replaced by C*) is parametrized by complex 
coordinates {z,w). We consider this curve in the type IIA background 

(IIA) X B X R2 X 5\ (2.7) 

and place a D4-brane wrapping x E and a D6-brane wrapping BxM.'^xS^. These branes 
intersect over E. Fermions on E are realized by massless modes of the 4-6 strings. The 
supersymmetry of the system ensures holomorphicity. The supersymmetries act trivially 
on the chiral fermions, which constitute a topological subsector of the complete system. 

Non-commutativity in this set-up is introduced by turning on a constant S-field along 
B, with holomorphic part 

B = \dz A dw. (2.8) 
A 

It is realized on the worldvolume of the D6-brane. 

By a chain of dualities presented in [5] this I-brane configuration relates to the back- 
ground 

(IIA) X X X S\ (2.9) 

with a D6-brane wrapping X x S^. This setup is appropriate for a computation of 
Donaldson-Thomas invariants DT{n, d), physically interpreted as BPS bound states of n G 
Hf){X,Z) = Z DO-branes and d e H2{X,'L) D2-branes to the D6-brane. The generating 
function of these invariants 

Zqu(i, A) = J] DT{n, d) e-^^e'^ * (2.10) 

n,d 
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Fig. 1: The I-brane configuration. A DJ^-hrane intersects with a D6- 
brane along a curve S. The 4-6 string degrees of freedom show up as 
free fermions on S. 



is closely related to the A-model topological string partition function on the toric manifold 
X with the complexified Kahler class t G H'^{X) 

Z,,,{t, A) = exp (^-^ - It . C2(X)^ Z^,{t, A). (2.11) 

Following the duality chain mentioned above, the (holomorphic) parameter A which 
initially specified a value of the S-field (2.8) acquires an interpretation of the topological 
string coupling constant A [5] . After summing over all bound states with p D4-branes as 
well, while weighting their contribution with a potential ^, the partition function of the 
final configuration reads 

Zi{UA)= J2 ef^^top(t + pA,A). (2.12) 

We identify this partition function with the I-brane partition function of the initial con- 
figuration (2.7). 

The above system can also be easily related to the supersymmetric gauge theories 
leading to a system of D4-branes spanned between NS5-branes. As shown by Witten [6], 
such a configuration engineers M = 2 supersymmetric gauge theories. We will elucidate 
this relation in much detail in section 5. 

B-field 
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The B-field quantizes the fermions on E. Let us first repeat the general arguments 
of [5]. The algebra A of open 6-6 strings on the D6-brane describes the interaction (as 
illustrated in figure 2) 

A^A^A, (2.13) 

and is explicitly non-commutative in the presence of a B-field. The S-field introduces a 
gauge field A on the D6-brane that couples to the open strings and quantizes the algebra 
of zero-modes of those strings [7, 8] . With a S-field given by 

B = \dz A dw (2.14) 
A 

the non-commutativity parameter is A. The complex coordinates z and w become non- 
commutative operators obeying 

[z,w]^X. (2.15) 

In case B — C x C we can identify this algebra with the Weyl algebra of differential 
operators 

A^Vc = {z,Xd,). (2.16) 

When we add to this system a D4-brane intersecting the D6-brane along the curve E, 
a 6-6 string acting on a 4-6 string can produce another 4-6 string (see figure 2) 

AxM^M. (2.17) 

This action endows the space of 4-6 open strings Ai with the structure of a module over 
the algebra A of 6-6 strings. Modules for the algebra of differential operators are called 
P-modules. 

To conclude, in the presence of a background flux, the chiral fermions on the 1-brane 
should no longer be regarded as sections of the spin bundle K^^'^. Instead they should be 
viewed as sections of a D-module. 

In the context of string theory it is worth stressing the range of parameters a' and A 
in which D-module description is valid. The string coupling A, which enters in the B-field 
flux as B = jdz A dw, plays an important role as quantization parameter. From the 
D-modulc point of view there seems to be no restriction on A, so one might hope that the 
P-module even captures non-perturbative information. However, in a particular system 
under consideration some restrictions on the values of A could arise that are related to the 
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radius of convergence of the partition function. Although we do make some additional 
remarks in Section 3 and in Section 6, we do not study these issues in this paper. 

On the other hand, the string scale ol does not play a fundamental role in the D- 
module. The P-module describes the topological sector of the intersecting brane configu- 
ration, which is realized in terms of massless modes of the I-brane system. Therefore the 
P-module description is valid only in the regime where ol is small (so that no massive 
modes interfere with our description). The most interesting case is of course when it is 
non-zero, as it provides a normalization factor for the worldsheet instanton contributions 
to the open 4-6 strings in the I-brane partition function (2.12). Section 5 clarifies this 
with an example. 




Fig. 2: The algebra A of functions on T, acts on the module M of free 
fermions. In the presence of a B -field the algebra A may be represented 
as a differential algebra, so that M becomes a D-module. 



V -modules and quantum curves 

We here introduce basic facts concerning P-modules and explain why they naturally 
describe I-branes. More details concerning theory of P-modules can be found in Ap- 
pendix B. 

"D-modules are defined as modules for the algebra of differential operators T>. In this 
paper we are interested in D-modules for the Weyl algebra V = (z, dz)- These are affine 
modules of rank 1 and can represented as 

V 



M 



V-P' 



where P is a linear differential operator P = aj(z)9* . The module Ad therefore captures 
solutions to the differential equation 

= 0, (2.18) 
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where takes values in some function space V, for example the algebra Oc of holomorphic 
functions on the complex plane C. "D-modules of rank 1 are cyclic, i.e. they are generated 
by a single element ^ E M., and so are of the form 

M = {D^ -.DeV}. (2.19) 

To be more precise, P-modules generated by the B-Held (2.8) depend on A and are 
known as 'D;^-modules [20] (when A is considered as a formal variable). Since all the 
differential modules we consider arc PA-modules, we often omit the subscript A. 

The D-module structure T> ■ M. C M gives a quantization of the semi-classical de- 
scription in equation (2.6). In particular, the rank 1 D-module 



is a quantization of the module 



>V=_-?i£^l (2.21) 

0{,,^} ■ F{z,w) ^ ' 



of functions on the curve defined by F{z, w) — 0. We therefore refer to the underlying 
differential equation P{z) = as a quantum curve. The I-brane set-up will obviously 
provide us with a rank 1 "D-module that represents a quantization of the I-brane curve 

Our notion of a quantum curve agrees with a notion of quantum spectral curves in the 
theory of (Hitchin) integrable systems. We discuss this relation shortly in Appendix C. 
Here we just give some examples of "D-modules and their interpretation in terms of quan- 
tum curves. 

Examples 

1) Take a linear partial differential operator on C, for example 

P = \zd^ - 1. (2.22) 

The differential equation = is solved by ^(z) = z^^^. So according to (2.19) the 
corresponding D-module can be represented as 

M = {z,\d,)z^'^. (2.23) 



^As a side remark notice that holonomic D-modulcs of dimension higher than 2 cannot be embedded 
in the 10 dimensions of string theory. Holonomic D-modules of dimension 2 are not related to the type 
II Calabi-Yau compactifications that we study in this paper, but could play a role in 4-dimensional 
Calabi-Yau compactifications in F-theory. 
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There are many equivalent ways of writing this module. For example, introducing 
^ — z^, the above differential equation is transformed into — with 

P = Xzd, - A - 1. (2.24) 

This new operator, as well as the solution to the new equation ^ = look different 

than before. Nonetheless, they represent the same P-module 

M = {z, Xd,) 2^+1/^ = {z, Xd,) z^l^. (2.25) 

This simple example illustrates how the formalism of D-modules allows to study solutions 
to partial differential equations independently of the way in which they arc written. 

An equivalent way to study "D-modules is in terms of fiat connections (see equation 
(B.IO)). The flat connection corresponding to P is given by 

Va = Xd^z - -dz, (2.26) 

z 

and determines ^{z) as a local flat section. It is a A-deformation of the degree 1 spectral 
cover 

E : w = -, (2.27) 

z 

with z,w &C*, together with the (meromorphic) 1-form 

A = -dz. (2.28) 

z 

This example enters string theory as the deformed conifold geometry describing the 
c—1 string. We will come back to it in section 4. 

2) All the modules that we will study in this paper are over C or C*. It is important 
that they may be of any rank though. Let us therefore also give a rank two example on 
the complex plane C. The second order differential equation 

= {X^d^ - z)^ (2.29) 

can be written equivalently as a rank two differential system 

P,,V', = 0, with P,, = Xd,)-{1 o)- (2-^°^ 
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Fig. 3: A second order differential operator P in Xdz defines a rank 2 
X-connection Va- The determinant 0/ Vq determines a degree 2 cover 
over C which is called the spectral curve S. 



Holomorphic solutions of this linear system are captured by the map 

that sends the two generators [(1, 0)*] and [(0, 1)*] to two independent (2- vector) solutions 
of P"^ = 0. The corresponding fiat connection 

VA^d,dz-j(^l l^dz (2.32) 

is a A-deformation of the degree 2 spectral cover (illustrated in figure 3) 

E : w^^ z, (2.33) 

with meromorphic 1-form rj = wdz\j:. Note that this one- form pushes forward to the 
connection 1-form, or Higgs field, 

II)". (2.34) 

in the basis {dz, wdz} of ramification 1-forms on the 2;-plane. Indeed, local sections of L 
push forward to local sections generated by 1 and w on the z-plane. Now, wdz ■ 1 = wdz 
and wdz ■ w — w^dz — zdz on E. 

We will discuss the string theory interpretation of this "D-module in detail in section 3. 
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2.3 Prescription 

Locally P-module on a general curve is just a system of linear difTerential equations, 
that changes from patch to patch. It is therefore natural to try to associate quantum 
invariants locally, and glue them with the help of the P-module transformations (which 
are a quantization of the usual coordinate transformations) . This is exactly the strategy 
we follow in the examples treated in the following sections. 

The simplest examples are deformations of affine curves of genus with a single marked 
point at infinity, i.e. consisting of a single patch near infinity. In these examples we 
associate a quantum invariant to the curve in the following way: 

We start with a differential equation P{z) — representing the quantum curve, where 
is the coordinate near infinity. Solutions to the differential equation P{z)'^{z) = 
form a module M., which is in particular O- module. We call this C- module W, and 
expect it to yield a subspace 

W C C{{z-^)). 

In that case we can, analogously as in the semi-classical case, turn W into a fermionic 
state and compute a tau-function. 

Since the l-brane configuration provides a "DA-module (in contrast to a "D-module) 
the resulting I-brane fermionic state |VV) is a A-deformation as well. We conjecture that 
its determinant computes the all-genus topological string partition function, when the 
appropiate quantum curve is chosen. 

In next sections we also discuss examples where we need to glue two local patches. 
Since each local patch yields a subspace of functions in the local variable, it is clear how the 
glueing should work: we need to insert a Fourier-like operator that relates the D-module 
on the first patch to the one on the second patch. As a result, the partition function turns 
into a correlation function which contracts the two corresponding fermionic states, with 
the insertion of the corresponding Fourier-like operator. (This is similar as in [4].) 

The above recipe doesn't tell us how to quantize a classical curve in a specific physical 
set-up. This is a very hard question in general. Moreover, it is not clear that the above 
prescription is independent of the chosen covering by local patches. Rather then provid- 
ing a general theory, in the following sections we analyze several important examples of 
spectral curves in string theory, and determine PA-modules that underlie their partition 
functions. Before we start with these examples though, let us explain the local procedure 
in more detail in two simple cases. 

Examples 
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1) Let's first explain the rank 1 case, with a "DA-module on C specified by the (mero- 
morphic) connection Va — — that may be triviahzed as 



When gx{z) is a holomorphic function on C that equals gx{0) — 1 at ^ = - in the notation 
of equation (A. 20) this is an element of r+ - this represents a pure gauge transformation 
on the disk, so that Va corresponds to a regular flat connection on C. 
For any gx{z) a fermionic section ^{z) of £ K^^'^ may be written as 



where ^(z) is a section of £ Cg) K^^^ with trivial connection d^- Flat sections "^{z) are 
defined by the differential equation 



They define a local trivialization of the bundle C with connection Va, and we will use 
them to translate the geometric configuration into a quantum state. 

A flat section for the trivial connection is given hj S^z) = 1. Wc associate the 
pair (C, dz) to the ground state \0) = z^ A z^ A z"^ A . . .. The gauge transformation gx{z) 
maps the trivial solution S(^) = 1 to ^^(2;) = gxi^), which transforms the vacuum into 
the fermionic state 



In other words, we build the quantum state by acting with the P-module generator 
^(2;) = gx{z) on the vacuum 



The state \W) is just the second quantization of the PA-module W. This state is non- 
trivial only when gx{z) is not a pure gauge transformation (which would correspond to a 
Krichever solution). In this situation the flat section diverges near z — 0, corresponding 
to a distorted geometry in this region. 

2) A degree n spectral curve E is quantized as a A-connection of rank n. This is 
equivalent to a "DA-module A4 that is generated by a single degree n differential operator 
P. As an Oc-module, though, M is generated by an n-tuple 



"^A^d.-gxiz) {d^gx^z)). 




\W)^gx\^)=gx{z)z'' ^gx{z)z' rK- 



W^Vx--^{z). 



(*(^),a,*(^),...,ar^*(^)). 
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where "^{z) is a solution of the differential equation = 0. In other words, this blends 
an n-vector of solutions to the linear differential system that the A-connection defines. 
We will name this Cc-module 

W = Oc-{'^iz),d,<i/iz),...,d:-'^{z)) C C((^-^)) 

(of course it contains the same elements as M.). This is the subspace we want to second 
quantize into a fermionic state |W). 

P has n independent solutions that differ in their behaviour at infinity. These 
solutions have an asymptotic expansion around z — oo that contains a WKB-piece plus 
an asymptotic expansion in A, and should thus be interpreted as perturbative solutions 
that hve on the spectral cover. We suggest that the asymptotic expansion of any solution 
can be turned into a fermionic state that captures the all-genus I-brane partition function. 
This partition function thus depends on the choice of boundary conditions near z^. 

Some of the WKB- factors will be exponentially suppressed near Zoo, while others grow 
exponentially. This depends on the specific region in this neighbourhood. The hues that 
characterize the changing behaviour of the solutions are called Stokes and anti-Stokes 
rays. Boundary conditions at infinity specify the solution up to a Stokes matrix: a solution 
that decays in that region can be added at no cost. 

This implies that the perturbative fermionic state we assign to a "D-module depends on 
the choice of boundary conditions. On the other hand, the P-module itself is independent 
of any of these choices and thus in some sense contains non-perturbative information and 
goes beyond the all-genus I-brane partition function. This agrees with the discussion in 
[27] . Nonetheless, the focus in this paper is on the perturbative information a "D-module 
provides. 

3 Matrix model geometries 

Hermitian one- matrix models with potential W{x) = Yl'j=o '^j^'' defined through 
the matrix integral 

In the large N hmit the distribution of the eigenvalues Aj of M on the real axis becomes 
continuous and defines a hyperelliptic curve. This curve is called the spectral curve of the 
matrix model. 
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In the 't Hooft limit this matrix model has a dual description as the B-model topo- 
logical string on Calabi-Yau geometries of the form [23] 



uv 



+ y^-W'{xy + f{x) = 0, (3.2) 



where f{x) — A/i Yl'j=o ^j^^ ^ polynomial in x of degree d—1. The hyperelliptic curve E 
modeling the local threefold equals the matrix model spectral curve, with 



N 



1=1 

The potential W{x) determines the positions of the cuts, containing the non-normalizable 
moduli, while the size of the cuts is determined by the polynomial /(x), comprising the 
normalizable moduli h^, . . . , 6^-2 and the log-normalizable modulus ha-i- 

This duality may be generalized by starting with multi-matrix models, whose spectral 
curve is a generic (in contrast to hyperelliptic) curve in the variables x and y. 

The I-brane picture suggests that the full B-model partition function on these Calabi- 
Yau geometries can be understood in terms of "D-modules. Even better, we will find that 
finite N matrix models are determined by an underlying D-module structure. 

In the past, as well as recently, these matrix models have been studied in great detail 
in several contexts. Most importantly for us, it has been realized that a central role is 
played by the string or Douglas equation 

[P,g] = A. (3.4) 

Here, P and Q are operators that implement multiplication with and differentiation with 
respect to a spectral coordinate. In a double scaling limit P and Q turn into differential 
operators. Physically, these critical models are known to describe minimal string theory. 

Already in [24, 25] an attempt has been made to understand this string equation in 
terms of a quantum curve in terms of the expansion in the parameter A. In Moore's 
approach this surface seemed to emerge from an interpretation of the string equation as 
isomonodromy equations. 

In topological as well as minimal string theory a dominant role is played by holomor- 
phic branes: either topological B-branes [4] or FZZT branes [26, 27] . Their moduh space 
equals the spectral curve, whereas the branes themselves may be interpreted as fermions 
on the quantized spectral curve. In particular, for (p, 1) minimal models the so-called Lax 
operator P has been interpreted as the quantization of the spectral curve. In these string 
theories it is possible to compute correlation functions using a l^i+oo-algebra [4, 28, 29] 
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that implements complex symplectomorphisms of the complex plane B in quantum theory 
as Ward identities. 

These advances strongly hint at a fundamental appearance of "D-modules in the theory 
of matrix models. Indeed, this section unifies recent developments in matrix models in 
the framework of section 2. Firstly, after a self-contained introduction in double scaled 
models we uncover the "D-module underlying the (p, l)-models. In the second part of this 
section we shift our focus to general Hermitian multi-matrix models, and unravel their 
X>-module structure. 

3.1 Double scaled matrix models and the KdV hierarchy 

Our first goal is to find the X>-modules that explain the quantum structure of double 
scaled Hermitean matrix models. This double scaling limit is a large N limit in which 
one also fine-tunes the parameters to find the right critical behaviour of the multi-matrix 
model potential. Geometrically the double scaling limit zooms in on some branch points 
of the spectral curve that move close together. Spectral curves of double scaled matrix 
models are therefore of genus zero and parametrized as 

Sp,, : yP + x'i + ... = 0. (3.5) 

The one-matrix model only generates hyperelliptic spectral curves, whereas the two- 
matrix model includes all possible combinations of p and q. These double scaled multi- 
matrix models are known to describe non-critical (c < 1) bosonic string theory based on 
the (p, g) minimal model coupled to two-dimensional gravity [30, 31, 32, 33, 34]. This 
field is therefore known as minimal string theory. 

Zooming in on a single branch point yields the geometry 

corresponding to the {p, 1) topological minimal model. This model is strictly not a well- 
defined conformal field theory, but does make sense as 2d topological field theory. For 
p = 2 it is known as topological gravity [35, 36, 37, 14]. 

All (p, q) minimal models turn out to be governed by two differential operators 

P = (Xd^f + Up^2{x){Xd,T-^ + ... + Uo{x), (3.6) 

Q = {Xd^y + Vg^2{x){Xd,y-^ + ... + vo{x), (3.7) 
of degree p and q respectively, which obey the string (or Douglas) equation 

[P,Q] = X. (3.8) 
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P and Q depend on an infinite set of times t — {ti,t2,t3, . . .), which are closed string 
couphngs in minimal string theory, and evolve in these times as 

a|-P=[(P^/^)+,P], (3.9) 
X-l-Q = [{P^/P)^,Q], (3.10) 

The fractional powers of P define a basis of commuting Hamiltonians^. This integrable 
system is known as the p-th KdV hierarchy and the above evolution equations as the KdV 
flows. 

The differential operator Q is completely determined as a function of fractional powers 
of the Lax operator P and the times t 

E + (3-13) 

j>i ^ 

j ^ mod p 

This implies that when we turn off all the KdV times except for ti = x and fix t^+i to be 
constant we find Q = Xd^. This defines the (p, l)-models 

P={\d,y-x, Q^Xd,. (3.14) 

One can reach any other (p, q) model by flowing in the times t. 

The partition function of the p-th KdV hierarchy is a tau- function as in equation (1.1). 
The associated subspace W E Gr may be found by studying the eigenfunctions 'ip{t, z) of 
the Lax operator P 

Pij{t,z) = zP'ip{t,z). (3.15) 

The so-called Baker function ip{t, z) represents the fermionic field that sweeps out the 
subspace W in the times t. 



^Notice that L = P^/p is a pscudo-diffcrcntial operator, having an expansion 

L = \d^ + lo{x) + h{x){Xd^)-^ + l2{x){Xd^)-^ + ..., (3.11) 
in negative powers of Xdx- This extended notion of a derivative is defined by the Leibnitz rule 



oo 

' n 



fe=o 



for any n G Z with (^) = n ■ . . . ■ {n — k + 1) / k\. It gives the derivatives with n < an interpretation of 
partial integration. L+ is the notation for the restriction to the positive powers of L. 
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If we restrict to the (p, l)-models the Baker function ip{x, z) can be expanded in a 
Taylor series 



Since iIj{x,z) is an element of W for all times, this defines a basis {vk{z)}k>o of the 
subspace W. In fact, it is not hard to see that the (p, 1) Baker function is given by the 
generalized Airy function 

i>{x,z) = e'(^^VzP-^ dw e ap/p+i ^ p+i , (3.17) 
which is normalized such that its Taylor components Vk{z) can be expanded as 

Vk{z)^z'{l + 0{X/z''+')) (3.18) 
The (p, 1) model thus determines the fermionic state 

\W) ^VoAviAv2A..., (3.19) 

where the Vk{z) can be written explicitly in terms of Airy-like integrals (see [14] for a nice 
review). The invariance under 

zP-WcW (3.20) 

characterizes this state as coming from a p-th KdV hierarchy. In the other direction, the 
state |W) determines the Baker function (and thus the Lax operator) as the one-point 
function 

^l^{t,z) = {t\^l>{z)\W). (3.21) 

In the dispersionless limit A ^ the derivative Xdx is replaced by a variable d, and 
the Dirac commutators by Poisson brackets in x and d. The leading order contribution 
to the string equation is given by the Poisson bracket 

{Po,go} = l, (3.22) 

where Pq and Qo equal P and Q eA X — 0. The solution to this equation is 

Po{d;t)=x (3.23) 
Qo{d;t) ^ y{x;t) (3.24) 
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and recovers the genus zero spectral curve Ep_g of the double scaled matrix model, 
parametrized by d. The KdV flows deform this surface in such a way that its singu- 
larities are preserved. (See the appendix of [27] for a detailed discussion.) 

Note that T,p^q is not a spectral curve for the Krichever map. The Krichever curve is 
instead found as the space of simultaneous eigenvalues of the differential operators 

[P,Q]=0, (3.25) 

that is preserved by the KdV flow as a straight-line flow along its Jacobian. In fact, there 
is no such Krichever spectral curve corresponding to the doubled scaled matrix model 
solutions. 

Wrapping an I-brane around Ej, ^ quantizes the semi-classical fermions on the spectral 
curve Ep q. The only point at infinity on E^ ^ is given by x — > oo. The KdV tau-function 
should thus be the fermionic determinant of the quantum state \W) that corresponds to 
this "D- module. In the next subsection we write down the "D-module describing the (p, 1) 
model and show precisely how this reproduces the tau-function using the prescription 
outlined in section 2. 

3.2 D- module for topological gravity 

We arc ready to reconstruct the "D-module that yields the fermionic state |yV) in equa- 
tion (3.19). For simplicity we study the (2, l)-model, associated to an I-brane wrapping 
the curve 

^(2,1)^ — X (with x, y e C). (3.26) 

Notice that this is an 2 : 1 cover over the x-plane. It contains just one asymptotic region, 
where a; ^ oo. Fermions on this cover will therefore sweep out a subspace W in the 
Hilbert space 

W C U{S') = C((y-i)), (3.27) 

the space of formal Laurent series in y~^. The fermionic vacuum |0) C T-C{S^) corresponds 
to the subspace 

\0)=y^/^Ay^/^Ay'>/^A..., (3.28) 

which encodes the algebra of functions on the disk parametrized by y and with boundary at 
y — oo. Exponentials in y~^ represent non-trivial behaviour near the origin and therefore 
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act non-trivially on the vacuum state. In contrast, exponentials in y are holomorphic on 
the disk and thus act trivially on the vacuum. 

The S-field B — jdxAdy quantizes the algebra of functions on into the differential 
algebra 

V^ = {x,Xd,). (3.29) 

Furthermore, it introduces a meromorphic connection 1-form A = j^ydx on S(2,i), which 
pushes forward to the rank two A-connection 

Va = A9,-(^° J) (3.30) 

on the base, parametrized by x. We claim that the corresponding I?;^-module M., gener- 
ated by 

P={\d,f-x, (3.31) 

describes the (2, 1) model. Let us verify this. 

Trivializing the A-connection Va in (3.30) implies finding a rank two matrix g{x) such 
that 

Wa = Xd^-g'{x)og-^{x). 

The columns of g define a basis of solutions ^{x) to the differential equation VA^ix) = 0. 
They are meromorphic fiat sections for Va that determine a trivialization of the bundle 
near x — oo. As the connection Va is pushed forward from the cover, ^(x) is of the form 



*(x) 



ijj{x) 
ip'ix) 



Independent solutions have different asymptotics in the semi-classical regime where 
X — > oo. In the (2, l)-model the two independent solutions ip±{x) solve the differential 
equation 

PM^) = {{Xd, f - x)^±{x) = 0. (3.32) 

Hence these are the functions ip+{x) — Ai{x) and ■0_(x) = Bi(x), that correspond semi- 
classically to the two saddles 



w± — ±^/x/\ 



1/3 
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of the Airy integral 

^(x) = dw e W^"^ 3 . (3.33) 

2'Ki J 

The 2^-module Ai can be quantized into a fermionic state for any choice of boundary 
conditions. Depending on this choice we find an (9(a;)-modulc W± spanned by hnear 
combinations of ip±{x) and of ip'^lx). The fermionic state is generated by asymptotic 
expansions in the parameter A of these elements. 

The saddle-point approximation around the saddle w± = ±^/x/\^^'^ yields 

V^±(x) ^ y-'/' e^'i + T.^^nX-{±y)-''^ 

To see the last step just recall the definition of vq{z) as being equal to the Baker function 
^^{x, z) evaluated &X, x — 0.^ A similar expansion can be made for with the result 

ip'±{x) - y^'^ e^* vi{±y). 

Note that both expansions in A are functions in the coordinate y on the cover. They 
contain a classical term (the exponential in 1/A), a 1-loop piece and a quantum expansion 
in \y^^. When we restrict to the saddle w = ^/x/X^^^, these series blend the into the 
fermionic state 

= My) A i^+{y) a y^My) a y^^y) a.... (3.34) 

Does this agree with the well-known result (3.19)? 

First of all, notice that the basis vectors x'^iplx) and x^^j\x), with /c > 0, contain in 
their expansions the function Vk{y) plus a sum of lower order terms in Vi{y) (with I < k). 
The wedge product obviously eliminates all these lower order terms. Secondly, the extra 
factor |/~^/^ factors just reminds us that we have written down a fermionic state. 

Furthermore, the WKB exponentials are exponentials in y and thus elements of the 
subgroup r+ of holomorphic functions that extend over a disc centered around y = Q, 
whereas the expansions Vk{y)/y^ are part of the subgroup r_ of functions that extend 
QYQi a disc centered at y — oc . (The definition of these subgroups and their action on the 

^Remark that x and appear equivalently in ■tp{x,z), while tp{x) and tl){x,z) only differ in the 
normalization term in z. 
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infinite Grassmannian can be found in Appendix A.) Up to normal ordening ambiguities 
this shows that the WKB part gives a trivial contribution to the fermionic state |W+). 
In fact, the tau- function even cancels these ambiguities. 
This shows that 



is indeed the same as in (3.19), when we change variables from 2; to y in that equation. 
Of course, this doesn't change the tau-function. 

So our conclusion is that the P-module underlying topological gravity is the canonical 
"D-module 



This "D-module gives the definition of the quantum curve corresponding to the (2, 1) model 
and defines its quantum partition function in an expansion around A = 0. Exactly the 
same reasoning holds for the (p, l)-model, where we find a canonical rank p connection on 
the base. It would be nice to be able to write down a P-module for general {p, g)-models 
as well. 

3.3 D-module for Hermitean matrix models 

"D-modules continue to play an important role in any Hermitean matrix model. In this 
subsection we are guided by [38] and [39, 40] of Bertola, Eynard and Harnad. 

We first summarize how the partition function for a 1-matrix model defines a tau- 
function for the KP hierarchy. As we saw before, such a tau-function corresponds to a 
fermionic state |W), whose basis elements we will write down. Following [38] we discover 
a rank two differential structure in this basis, whose determinant reduces to the spectral 
curve in the semi-classical limit. This "D-module structure is somewhat more complicated 
then the "D-module we just found describing double scaled matrix models. 

We continue with 2-matrix models, based on [40] . Instead of one differential equation, 
these models determine a group of four differential equations, that characterize the V- 
module in the local coordinates z and w at infinity. The matrix model partition function 
may of course be computed in either frame. 

1-matrix model 




(3.35) 



M = 



(3.36) 
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Let us start with the 1-matrix model partition function (3.1). By diagonahzing the 
matrix M the matrix integral may be reduced to an integral over the eigenvalues Aj 

Zn^ f n^A, A(A)2e-xE.H.(A.), (3.37) 

with the Vandermonde determinant A(A) = nj<j(-^i ~ '^i) — det(A^~^). The method of 
orthogonal polynomials solves this integral by introducing an infinite set of polynomials 
Pk{x), defined by the properties 

Pk{x)=x\l + Oix-')), (3.38) 
J dx pk{x) pi{x) e-i^(-) = hkS,,i. (3.39) 

The normalization of their leading term determines the coefficients /i„ G C. Since the Van- 
dermonde determinant A{x) is not sensitive to exchanging its entries xj ^ for pj^i{xi), 
substituting A{x) = det(pj_i(a;j)) turns the partition function into a product of coeffi- 
cients 

iV-l 

Zj, = N\l[ hk. (3.40) 

fc=0 

With the help of orthogonal polynomials the large N behaviour of may be studied, 
while keeping track of corrections. 

In this discussion the orthogonal polynomials arc relevant since they build up a basis 
for the fermionic KP state. In an appendix of [38] it is shown that one should start at 
t = with a state | Wq) generated by the polynomials Pk{x) for k > N 

I Wo) = PNix) A Pn+i{x) a Pn+2{x) a.... (3.41) 

Notice that the vector Pn{x) thus corresponds to the Fermi level and defines the Baker 
function in the double scaling limit. Acting on them with the commuting flow generated 

by 

r+ = |(?(t) =e^">i^*"^"} (3.42) 

defines a state \ Wt) — \g{t)Wo) at time t, which allows to compute a tau-function at time 
t. If the coefficients Uj in the potential W(x) are taken to be Uj — Uj^^ +tj, this r-function 
equals the ratio of the matrix model partition function Zn at time t divided by that at 
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Multiplying the orthogonal polyonomials by exp{—-^W{x)) doesn't change the 
fermionic state W — Wo in a relevant way, since this factor is an element of r+. To 
find the right "D-module structure, it is necessary to proceed with the quasi-polynomials 

1 



--Pke 



which form an orthonormal basis with respect to the bilinear form 



(3.43) 



(3.44) 



It is possible to express both multiplication by x and differentiation with respect to x 
in terms of the basis of V'm's- The Weyl algebra {x, Xd^) acts on these (quasi)-polynomials 
by two matrices Q and P 



xipk{x) = ^Qkiipi 

1=0 

oo 

Xdxi^kix) = '^Pki'ipiix), 



(3.45) 
(3.46) 



1=0 



and the space of quasi-polynomials ipk is thus a "DA-module. 



Notice that we anticipate that the D-module possesses a rank two structure, since 
we started with a flat connection A = jydx on an I-brane wrapped on a hyper-elliptic 
curve. Now, the matrices Q and P only contain non-zero entries in a finite band around 
the diagonal. The action of dx on the semi-infinite set of iljk{x)^s can therefore indeed be 
summarized in a rank two differential system ([38] and references therein) 



An{x) 



(3.47) 



where A]\f{x) is a rather complicated 2 x 2-matrix involving the derivative W of the 
potential and the infinite matrix Q: 



An{x) = ^W'{x) 



"-10" 




1 


+ 7n 



-\r{Q,x)N,N-l W;iQ,x)N,N 

N-1,N-1 N-1,N 



(3.48) 



with 



W'(Q,x)^ 



W'{Q) - W'{x) 
Q - X 



and 



h 



N 



h 



N-1 



(3.49) 
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Equation (3.47) is thus the rank two A-connection defining the "DA-module structure on 
W that we were searching for! As a check, the determinant of this connection reduces to 
the spectral curve in the semiclassical, or dispersionless, hmit [38]: 



N-l 



= det (yl2x2 



(3.50) 
(3.51) 



jj 



(To make the coefficients in the above equation agree with (3.2), we rescaled y 1//2.) 
In conclusion we found the "D-module structure underlying Hermitean 1-matrix models. 

Remark that in the N ^ oo limit we expect that the hyperelliptic curve defining the 
B-model Calabi-Yau (3.2) emerges from Eat. Indeed, in the 't Hooft limit Q corresponds 
classically to the coordinate x on the curve, whereas quantum-mechanically it is an oper- 
ator whose spectrum is described by the eigenvalues Aj of the infinite matrix M. In the 
large limit we can therefore replace the matrix Qij in the definition for Eat by XiSij. 

Note as well that we can rewrite the rank two connection for the vector (-^at, i/^'nY ^ 



A9x 







" 1 " 






V'jv(^) 




-det(^jv(a;)) + Ar XZ 







(3.52) 



at least when tT:{AN{x)) — 0, with Y and Z some derivatives of entries of An{x). This 
brings the A-connection in the familiar form of section 2. In the next subsection we clarify 
the differential structure in a simple example. 

2-matrix model 

Let us first say a few words on the X>-module structure underlying multi-matrix models, 
which capture spectral curves of any degree in x and y [39, 40]. The partition function 
for a two-matrix model, with two rank N matrices Mi and M2, is 



-I 



DM1DM2 e 



j-Tr{Wi{Mi)+W2{M2)-MiM2) 



(3.53) 



where Wi and 14^2 s^ie two potentials of degree di + 1 and ^2 + 1. Choosing W2 to be 
Gaussian reduces the 2-matrix model to a 1-matrix model. The 2-matrix model is solved 
by introducing two sets of orthogonal polynomials Trk{x) and (Tk{y). Again it is convenient 
to turn them into quasi-polynomials 



(3.54) 



28 



obeying the orthogonality relations 

I dxdy ipk{x)4>i{y)(^^ = hkhi- (3.55) 

Multiplying with or taking a derivative with respect to either x oi y yields (just) two 
operators Q and P (and their transposes because of (3.55)), that form a representation 
of string equation [P, Q] = 0. Since Q is only non-zero in a band around the diagonal of 
size d2 + l and P of size di + 1, the quasi-polynomials may be folded into the vectors 

— * , — * , 

1p^[lpN,---,1pN-d2] , = [0JV, . . . >JV-<ii] . (3.56) 

Any other quasi-polynomial can be expressed as a sum of entrees of these vectors, with 
coefficients in the polynomials in x and y. These vectors arc called windows. The differ- 
ential operators \dx and \dy respect them, so that their action is summarized in a rank 
d2 + 1 resp. rank di + 1 A-connection 

XdJix) = A,{x)if{x), \dy${y) = A2{y)${x). (3.57) 

This we interpret as two representations of the 'D\-m.odv\c underlying 2-matrix models. 
Indeed, [39] proves that the determinant of both differential systems equals the same 
spectral curve S, in the limit A — when we replace \dx y and \dy x. The defining 
equation of E is of degree rfi + 1 in x and of degree d2 + 1 in y. 

In fact, it is useful to introduce two more semi-infinite sets of quasi-polynomials i^J^y) 
and (pj^ix), as the Fourier transforms of ipk{x) and (pkiv) respectively. The action of the 
Weyl algebra on them may be encoded as the transpose of the above linear systems. The 
full system can therefore be summarized by (compare to (4.7)) 

X-axis : {i/Jkix), ^^(x)}, Va = Xd^ - Ai{x), (3.58) 
l/-axis: {0fe(|/), 3/^^(1/)}, V x = Xdy - A2{y) . 

Moreover, the matrix model partition function can be rewritten as a fermionic correlator 
in either local coordinate 

Zn^ f UdXldX^ A{X')A{X') e-i^>(^^(^^)+^^(^?)-^«^? (3.59) 

N-1 N-1 
k=0 k=0 

with respect to the bilinear form in (3.44). 
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Furthermore, Bertola, Eynard and Harnad study the dependence on the parameters 
u'p and u^p appearing in the potentials Wi and W2- Deformations in these parameters 

— * — * 

leave the two sets of quasi-polynomials invariant as well. On and they act as matrices 
Uf^ and Uf \ This yields the 2-Toda system 

d^i.)Q = -[g, Uf] d^wP - -[P, f^j'^] (3-60) 
9^(.)g = [Q, C/f ] P = [P, (3.61) 

3 3 

In [39] it is proved that the linear differential systems (3.57) arc compatible with these 
deformations, so that the parameters u~p and u~p in fact generate isomonodromic de- 
formations. This shows precisely how the non-normalizable parameters in the potential 
respect the central role of the Vx-modvle (3.57) in the 2-matrix model. 

3.4 Gaussian example 

Let us consider the Gaussian matrix model with quadratic potential 

W{x) = ^, (3.62) 

that is associated to the spectral curve 

y2 ^ - 4//^ (3.63) 

in the large N limit. In the Dijkgraaf-Vafa correspondence this matrix model is thus dual 
to the topological B-model on the deformed conifold geometry. 
The Hermite functions 

ijlix) = -^e-^^Hl{x), with 



X 




form an orthogonal basis for this model. Their inner product is given by 



/ 



The partition function of the Gaussian matrix model can be computed as a product of 
the normalization constants hk. Using the asymptotic expansion of the Barnes function 
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G2{z), that is defined by G2{z + 1) = r{z)G2{z), the free energy can be expanded in 
powers of A 

N-l 

(3.64) 



^iv = logn^fc = log(<^2(A^ + l) 

k=l 
1 //i\2 



(27r)^/2 



2<7 



5=2 



2^(2^-2) V/^y 



where are the BernouUi numbers and ^ — NX. 

The derivatives of the Hermite functions are related as 



A 



d_ 
dx 







-x/2 






-Vkx 



kX 



^l{x) 



So, according to the previous discussion, the I>A-module connection is given by 

d 



^ d - / s . d 



dx 



x/2 -VPfX 
^/WX -x/2 



(3.65) 



(3.66) 



Here we choose ip = [ipN, i'N-iY as window. In the large N limit the determinant of this 
rank two differential system indeed yields the spectral curve (3.63) with /i — XN. 

Instead of using -0^ and 'ipk-i as a basis, we can also write down the differential system 
for and its derivative ip'ki^) — Since this derivative is a linear combination 

of V'fe-i and xipk{x) (as we saw above), it is equivalent to use this basis to generate the 
fermionic state W. We compute that 



A^ 

dx 



X 





AA^ 



A/2 



1 





.X 



(3.67) 



The spectral curve in the large N limit hasn't changed. Notice that in this form it is 
clear that the rank 2 connection is the push-forward of the connection A — \ydx on the 
spectral curve — x^ — 4// to the C-plane, up to some A-corrections. 

In the double scaling limit the limits N ^ oo and A — > are not independent as in the 
't Hooft limit, but correlated, such that the higher genus contributions to the partition 
function are taken into account. In terms of the Gaussian spectral curve this limit implies 
that one zooms in onto one of the endpoints of the cuts. The orthogonal function ipj^ix) 
turns into the Baker function ^p{x) of the double scaled state W. 

In the Gaussian matrix model this is implemented by letting x — > y^ + ex, where e is 
a small parameter. So the double scaled spectral curve reads 



(3.68) 
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while the differential system reduces to 







"01" 




ipix) 


^|;'{x) 




X 




^|;'{x) 



(3.69) 



This is indeed the P-module corresponding to the (2, l)-model. 

4 Conifold and c = 1 string 

The free energy (3.64) of the Gaussian matrix model pops up in the theory of bosonic 
c = 1 strings. This c = 1 string theory is formulated in terms of a single bosonic coordinate 
X, that is compactified on a circle of radius r in the Euclidean theory. A critical bosonic 
string theory (with c = 26) is obtained by coupling the above CFT to a Liouville field (p. 
The Liouville field corresponds to the non-decoupled conformal mode of the worldsheet 
metric. The local worldsheet action reads 

^ / (^(^^)' + (^'^)' + ^^^^ + ^'^^) ' 

where the coupling is seen as the worldsheet cosmological constant. In the Euclidean 
model there are only two sets of operators, that describe the winding and momenta modes 
of the field X. These vertex and vortex operators can be added to the action as marginal 
deformations with coefficients tn and 

Just like in c < 1 minimal string theories (the (p, g)-models of last section), the parti- 
tion function of the c — 1 string is first computed using a dual matrix model description 
[41]. At the self-dual radius r = 1 it agrees with the Gaussian matrix model partition 
function in equation (3.64), where A now plays the role of the c — 1 string coupling 
constant. 

The matrix model dual to the c — 1 string is called matrix quantum mechanics. 
This duahty is reviewed in much detail in e.g. [42, 43, 44]. Matrix quantum mechanics 
is described by a gauge field A and a scalar field M that are both N x N Hermitean 
matrices. The momentum modes of the c — 1 string correspond to excitations of M, 
whereas the winding modes are excitations of A. If we focus on the momentum modes, 
the (double scaled) matrix model is governed by the Hamiltonian 

2 V J 
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Let us focus on solutions that depend purely on the eigenvalues Aj of M. The Hamil- 
tonian may be rewritten in terms of the eigenvalues as 

H-lA-\X)Y,(-X'^,-Xi) A(A), 

where A (A) is Vandermondc determinant. It is convenient to absorb the factor A in the 
wavefunction solutions, making them anti-symmetric. Hence, the singlet sector of matrix 
quantum mechanics describes a system of N free fermions in an upside-down Gaussian 
potential. 

To describe the partition function of the c = 1 model it is convenient to move over to 
light-cone coordinates \± = \ ± p, so that elementary excitations of the c = 1 model are 
represented as collective excitations of free fermions near the Fermi level 

A+A_ = Pi. (4.1) 

When we restrict to A± > 0, scattering amplitudes can be computed by preparing asymp- 
totic free fermionic states {t\ and \t) at the regions where one of A± becomes very large. 

In this picture the generating function of scattering amplitudes has a particularly 
simple form. It can be formulated as a fermionic correlator [45] 

Z^{t\S\t), (4.2) 

where the fermionic scattering matrix S G GL[oo,C) was first computed in [46]. More- 
over, in [47] (see also Chapter V of [44]) and later in [4] it is noticed that S just equals 
the Fourier transformation 

(5V')(A_) = J dA+ei^-^+V(A+). (4.3) 

In the next section we show that this follows naturally from the perspective of D-modules. 

The result (4.2) shows that c — 1 string theory is an integrable system, just hke the 
(p, g)-models in the last section. Since it depends on two sets of times this integrable 
system is not a KP system. Instead, the above expression defines a tau function of a 
2- To da hierarchy. 

Notice that the Fermi level (4.1) is a real cycle on the complex curve 

E : zw = ji, (4.4) 

which is a different parametrization of the spectral curve — — /i oi the Gaus- 
sian 1-matrix model. In the revival of this subject a few years ago, a number of other 
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matrix model interpretations have been found. This includes a duality with the Her- 
mitean 2-matrix model, which makes the 2-Toda structure manifest [48], a Kontsevich- 
type model [49, 50] at the self-dual radius, and a so-called normal matrix model [51, 52], 
that parametrizes the dual real cycle on the complex curve E. Let us also mention that 
the well-known duahty of the c—1 string with the topological B-model on the deformed 
conifold [53] , that follows, with a detour, from the more general Dijkgraaf-Vafa correspon- 
dence. 

4.1 D- module description of the c—1 string 

This paragraph reproduces the c—1 partition function (4.2) from a I>-module point 
of view. The discussion continues the line of thought in section 5.5 of [4] and in [5]. 

As we have just seen, the c—1 string is geometrically characterized by the presence 
of a holomorphic curve in C x C defined by 

Ec=i : zw — II. 

Let us consider an I-brane wrapping the curve Ec=i. When we assume z as local coordinate 
the curve quantizes into the differential operator 

P = -Xzd, - II. (4.5) 

It is amusing that the differential operator P appears as a canonical example in the theory 
of "D-modules (see e.g. [11]) in the same way as the c—1 string is an elementary example 
of a string theory. 

We recognize this example from section 2, where a "D-module was associated to the 
differential operator P. However, now it is important not to forget that there are two 
asymptotic points 2^00 and Woo- Let us call their local neighbourhoods Uz and Uu,, as local 
coordinates are z and w respectively. At both asymptotic points the I-brane fermions 
will sweep out an asymptotic state. The quantum partition function should therefore be 
constructed from two quantum states. 

Before constructing these states for general A, let us first consider the semi-classical 
hmit A — > 0. In this limit the I-brane degrees of freedom are just conventional chiral 
fermions on Ec=i. The genus 1 part Ti of the free energy is obtained as the partition 
function of these semi-classical fermions. It can be computed by assigning the Dirac 
vacuum 

\0)z^z'/'Az'/'Az'/'a... 
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to Uz and likewise the conjugate state 

to Uw To compare these states, we need an operator 5" that relates z to 1/z. The semi- 
classical partition can then be computed as a fermionic correlator ^(0 15*10) 2, with the 
result that 



— w 



{OlS^z^Hl^'+'Z^ (4.6) 



fc>0 



Using ^-function regularization we find that this expression yields the familiar answer 
^1 = -nlog/x. 

In order to go beyond 1-loop, we should think in terms of X>-modules. Let us for 
a moment not represent their elements in terms of differential operators yet. In both 
asymptotic regions we then find the P-modules 

U^: M = V/VP, with P = zw - 

U^: M^V/VP, with P^wz-pi + X. 

Notice that the Weyl algebra V = {Zjw), with the relation [z,w] = A, acts on mono- 
mials 2;^ and w'^ in the module Ai as 

z{z'') = z''^^ z{w'') = (^Xd^ + w'' 

w{z'') = (-Aa, + ^) z'' w{w'') = w''+\ 

Here, we just used the relation VP = and wrote the elements in the basis {z^, \ k E 
Z} of A basis of a representation of Ai on which z and w just act by multiplication 
by z resp. differentiation with respect to z is given by 

v'^{z) = z'' ■ z-^"/^, 



Indeed, differentiation with respect to z clearly gives the same result as applying w. 
Moreover, multiplying by z gives 

z ■ vl{z) ^xjdw e"'"'/^^ = (// + A(A; - 
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Similarly, in the module M one can verify that 

\ w/ 
ziw"^) = (^Xd, + z"" z{z'') = z''+\ 

Hence in the representation of M defined by 

vl{w) ^ J dz e^"'/^ z'' ■ z-"/^, 

w and dy, act in the usual way. 

Since we moved over to representations of the P- module where the differential operator 
acts as we are used to, the 5" transformation, that connects the 11^ and the 11^ patch and 
thereby exchanges z and w, must be a Fourier transformation. This is clear from the 
expressions for the basis elements w and w: S interchanges v^{z) with Vk{'^)-> ^^'^ "^ki^) 
with Vk{w). In total we thus find the D-module elements 

U. : vl < (4.7) 
Uw '■ VLki V.k 

Representing the X>-module in terms of differential operators of course gives the same 
result. A fundamental solution of P'i{z) ~ is ^{z) = z'^^^^, so that acting with 
V = {z,dz) on '^{z) gives the elements vl'm. M.. Likewise, we reconstruct the elements 

from the fundamental solution of P^{w) — 0. Since V — {z, d^) and V = {w, dyj) are 
related by a Fourier transform, an element Vk of the I>-module in one asymptotic region is 
represented by its Fourier transform in the opposite region. This reproduces all elements 
in (4.7). 

A A-expansion of the "D- module element using for example the stationary phase 
approximation, yields as zeroth order contribution 



e 



jJi \ 



(-) 



while the subdominant contribution is given by 



2nXfj, 
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So in total we find that 



qu 



This summarizes the contributions that we found before: the genus zero w^^^ and genus 
one ii^^^/'^w~^~^ results, plus the higher order contributions that are collected in ■0, 



qu- 



The all-genus partition function Z of this I-brane system can be easily computed 
exactly. Schematically it equals the correlation function 



{y^u,\s,\w,), 



where the ^'-matrix implements the Fourier transform between the two asymptotic 
patches. Similar to the arguments in (the appendices of) [47] and [4]^ wc find that the 
result reproduces the perturbative expansion of the free energy as in equation (3.64). For 
completeness let us review the argument by comparing v^(w) with v^{w). 

Notice that 'if^.iw) almost equals the gamma-function r(^) = f^dt e~* t^~^. Indeed, 
let us take the integration contour from —ioo to ioo and choose the cut of the logarithm 
to run from to oo. Then 



^^(-) = ( ^ 



fc+i- 



fc+i- 



lOO 





dz' e'' — 
\ w 



dz' e 



iz' ^(fe-f)logz 



-oo 




dz' e''' e('=-A)'°g^ 



/■oo 

' + / dz' e''' 
Jo 

^0 



,(fc-^)logz' 



dz' e 



tz' g(A;-^)logz' 



where we moved the contour along the positive imaginary axis. A change of variables and 
using that log{iz' — e) — log 2;' — 3i7r/2 and log{iz' + e) — logz' + in/2, for e small and 
real, then yields 



w 



^^i(fc+l_£)/2 



-37ri(fc+l-^)/2 



r k + 1- 



which is the same as the theory of type II result in the appendix of [47]. Ignoring the 
exponential factor (which will only play a role non-perturbatively) , we find that the free 
energy T equals the sum 

T{X, ^) = J2{k + l-^)logX + logr{k + l-^). 



k>0 



^The argument presented in the appendix of [4] is not fully correct. The proper argument (as shown 
below) recovers a slightly different prefactor in front of the Gamma-function, related to the doubling in 
the appendix of [47]. 
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It obeys the recursion relation 

which is known to be fulfiUed by the c = 1 string (see for example Appendix A in [55]), 
up to a term — | log(27rA) that can be taken care of by normalizing the functions Vf.. The 
same result is found when analyzing the function Vk- 

This concludes our discussion of the c = 1 string. It is the first "D-module example 
where we see how to handle curves with two punctures. The physical interpretation of 
the 1-brane set-up furthermore provides a check of our formalism. Moreover, this example 
agrees with the claim that the D-module partition function should be invariant under 
different parametrizations. Both the representation as c = 1 curve, Sc=i : zw = fi, and 
that as a Gaussian matrix model spectral curve, Smm : 2/^ = a;^ + /x, yield the same 
partition function. 

5 Seiberg-Witten geometries 

Many times M = 2 supersymmetric gauge theories proved to provide an important 
theoretical framework for testing new ideas in physics. It should be fair to say that the 
most important advances in this context are the solution of Seiberg and Witten in terms of 
a family of hyperelliptic curves, as well as the explicit solution of Nekrasov and Okounkov 
in terms of two-dimensional partitions. In what follows we will provide a novel perspective 
on these results, by wrapping an I-brane around a Seiberg-Witten curve. The S-field on 
the I-brane quantizes the curve, and a fermionic state is obtained from the corresponding 
"D-module. As we will see, this state sums over all possible fermion fluxes through the 
Seiberg-Witten geometry, and may be interpreted as a sum over geometries. First we 
briefly review the Seiberg-Witten and Nekrasov- Okounkov approaches. 

The solution of the U{N) Seiberg-Witten theory is encoded in its partition function 
Z{ai, \, A), which is a function of the scale A, the coupling A and boundary conditions 
for the Higgs field denoted by for i = 1, . . . , iV (with ai — for the SU {N) theory). 
The partition function is related to the free energy as 

Z(a„A,A) =e^ = e^r=oA^-^^.K.A). (5.I) 

In the above expansion J^o is the prepotential which contains in particular an instanton 
expansion in powers of A^-^, while higher ^^'s encode gravitational corrections. The U (N) 
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Seiberg-Witten solution identifies the Oj's and the derivatives of the prepotential 
as the Ai and Bi periods of the meromorphic differential 

1 dt ^, 

on the hyperelliptic curve 

■ ^""{t + = PNiy) = ^{v - Ui). (5.3) 

i=l 

Despite great conceptual advantages, extracting the instanton expansion of the prepo- 
tential from this description is a non-trivial task. However, an explicit formula for the 
partition function, encoding not only the full prepotential but also entire expansion in 
higher Tg terms, was postulated by Nekrasov in [54]. Subsequently this formula was de- 
rived rigorously jointly by him and Okounkov in [55] and independently by Nakajima and 
Yoshioka in [56, 57]. For U{N) theory this partition function is given by a sum over N 
partitions R — {R{i), . . . ,R(n)) 

Z{a„X,A) = ZP-*(a„A)^A2^l^l/xJ(a„A), (5.4) 



R 



where 



Hs[ai,A)= \\ — r , (5.5j 

ZP-*(a,, A) = exp ( ^ 7,(a, - a,-. A)) . (5.6) 

The function 7a (a;. A) is related to the free energy of the topological string theory on the 
conifold, and its various representations and properties are discussed extensively in [55] 
in Appendix A. The vevs Oj are quantized in terms of A, so that for pi e Z, 

ai = \{Pi + Pi), Pi = — . (5.7) 

The approach of [54] is based on the locahzation technique in presence of the so-called 
Q-background. In general this background provides a two-parameter generalization of 
the prepotential: the coupling A is replaced by two geometric parameters ei and €2- The 
prepotential, as given above, is recovered for A = ei = —62. By the duality web in 
§?? supersymmetric gauge theories are related to intersecting brane configurations. The 
Nekrasov- Okounkov solution must therefore have an interpretation in terms of a quantum 
Seiberg-Witten curve, where A plays the role of the non-commutativity parameter. 
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5.1 Dual partition functions and fermionic correlators 

For a relation to the I-brane partition function (2.12), it is necessary to consider the 
dual of the partition function (5.4). This is introduced in [55] as the Legendre dual 

Z^(e,p,A,A)= Yl Z{\{p, + p,),\K)ei^^^^^K (5.8) 

An important observation of Nekrasov and Okounkov is that this dual partition func- 
tion can be elegantly written as a free fermion correlator. This is a consequence of the 
correspondence between fermionic states and two-dimensional partitions described in Ap- 
pendix A. For U{1) there is no difference between the partition function and its dual and 
both can be written as 

Z^(,)(p,A,A) = (p|e-^"^A2^°eX"-|p), (5.9) 

where \p) is the fermionic vacuum whose Fermi level is raised by p = a/ A units and Lq 
measures the energy of the state. A version of the boson-fermi correspondence implies 
the following decomposition 



R 



n — m J--*- h[lJ) 



in terms of partitions R, where fiR is the Plancherel measure 

Rm - Rn + n-m 

n — m ^ 

l<m<n<oo Dei? 

which can be written equivalently as a product over hook lengths /i(n). 
For general N the dual partition function (5.8) looks very similar 

^^(iv)(6;P,A,A) = (p|e-i°^e^^^A2^"ei"-|p), (5.12) 

however, now this expression is obtained by blending N free fermions ■0^'^ into a single 
fermion ■0, as explained in Appendix A. In particular H^. — ^ ^(r+1/2) mod Ar'0r'0-r) 
while the bosonic mode q;_i arises from the bosonization of the single blended fermion 
■0. In formula (5.10) the Plancherel measure of a blended partition R can be decomposed 
into N constituent partitions as 



/i^(ai,A), (5.13) 



with jjL^ and Z"^^^^ given in (5.5) and (5.6). When read in terms of the N twisted fermions 
the correlator (5.12) involves a sum over the individual fermion charges pj. 
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Our aim in this section is to derive the above fermionic expressions for the dual par- 
tition function from the perspective of this paper. In the next subsections we will see 
how canonically quantizing the Seiberg-Witten curve in terms of a D-module elegantly 
reproduces to the fermionic correlators (5.9) and (5.12). 

5.2 Fermionic correlators as P- modules 

In this section we compute the I-brane partition function for U{N) Seiberg-Witten 
geometries. We start with the simpler U{1) and U{2) examples and then generahze this 
to U{N). As a first principal step we notice that the U{N) Seiberg-Witten geometry 

N 

Hsw ■ h-'^it + r^) = Pn{v) = \[{v - Ui), (5.14) 

i=l 

can be rewritten as 

{Pm{v) - A''t){PN{v) - A^r^) = A^^. (5.15) 

This shows that the Seiberg-Witten surface may be seen as a transverse intersection of a 
left and a right half-geometry defined by 

Ei : A^t = Pn{v) resp. : A^r^ = Pn{v), (5.16) 

which are connected by a tube of size A^^. The left geometry parametrizes the asymptotic 
region where both t — > oo and v — > cxo, whereas the right geometry describes the region 
where v — > cxo while t ^ 0. This is illustrated in figure 4. 

Next we wish to associate a subspace in the Grassmannian to both half Seiberg-Witten 
geometries. This will be swept out by a fermion field on the curve that couples to the 
holomorphic part of the S-field 

B = \dsAdv (5.17) 
A 

Since this 5-field quantizes the coordinate v into the differential operator Xdg, any sub- 
space in this section is a P-module for the differential algebra 

Dc*^{t,Xds). (5.18) 

The free fermions on the Seiberg-Witten curves couple to the gauge field A — jrjsw- 
This determines their flux through the Ai cycles of the Seiberg-Witten geometry as 

Pi^ \ 1 Vsw- (5.19) 
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Fig. 4: The right-half Seiberg- Witten geometry is distorted around the 
asymptotic point {t ^ 0,v oo). A fermion field on the quantized 
curve can be described as an element of a V-module, and sweeps out a 
state |yV) at the -boundary where t — ^ oo. 



The flux leaking through inflnity is p = which is zero for SU{N). A fermion fleld 

with fermion flux p at inflnity, will sweep out a fermionic state in the pth Fock space. The 
parameters C,i = f]sw are dual to the fermion fluxes. Notice that in the perturbative 
regime Pi can be written as a A-expansion 

Xpi^Ui + 0{X). (5.20) 

Since both half Seiberg- Witten geometries are distorted near v — oo (see figure 4), 
while a fermionic subspace can be read off in the neighbourhood where v is finite, both 
half-geometries parametrize a subspace oi C{{v)): 

Si, S^ C C{{v)). (5.21) 

The trivial geometry corresponds to a disk with origin at v = oo, whereas its boundary 
encloses the point v — 0. The vacuum state is therefore given by 

|0) = A v-^/^ A y-^/^ A.... (5.22) 

Exponentials in act trivially (as pure gauge transformations in r+) on this state, 
whereas exponentials in v transform the vacuum into a non-trivial fermionic state. 
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Finally, the partition function is recovered by contracting the left and the right 
fermionic state. Note that s — — log i is a local spatial coordinate on both half Seiberg- 
Witten geometries, which tends to — oo on the left and to +00 on the right. This makes a 
huge difference with the c—1 geometry discussed in [4, 5] , where the local coordinate is 
the exponentiated coordinate, which on the left is the inverse of that on the right. While 
in that example a non-trivial S'-matrix is required to identify the left and right half- 
geometries, here we can just glue the fermionic states using the classic Hamiltonian Lq. 

Let us now find these quantum states! 

theory 

The U{1) Seiberg-Witten curve is embedded in C* x C as 

A{t + r^) ^v-u, {t^e' eC*, V eC) (5.23) 

where w G C is a normalizable mode. This geometry may be factorized into a left and a 
right geometry 

El: v^At + u and T,r : v^Ar^ + u, (5.24) 

that intersect transversely with degeneration parameter A^. 

The symplcctic form B = jds A dv quantizes both half geometries into ©A-modules 
on a punctured disc C^, parametrized by t. We claim that these are characterized by the 
U{1) A-connections 

Vl = -Ai9t + Ai + Ap and V r = Xtdt + Ar^ + Xp. (5.25) 

These are just the canonical quantizations of the classical Seiberg-Witten geometries, 
where additionally u is quantized into Xp, with p e Z. They yield the linear differential 
equations 

PL^l{t-p) = {-Xtdt + M + Xp)i,l{t-p) = Q, (5.26) 
Pni^iit-p) = {Xtdt + At-^ + Xp) ^l{t-^;p) = 0. (5.27) 

The "DA-modules are of the canonical form 



Ml/r = ^ , (5.28) 

'-'X ■ ^L/R 



and are generated by the solutions 

^l{t]p)^tPe^* and 'ijj^(t;p) ^ t-Pe^^~\ (5.29) 
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Prom the discussion in Appendix A it follows that the factor f'^ acts on the right 
Dirac vacuum by raising the Fermi level into \p), while the exponent of translates 
to the exponentiated q;_i operator. With an analogous statement for the left state, the 
modules M.l/r translate into the Bogoliubov states 

(WlI = (p|e^"i and | Wi?) = e^"-i (5.30) 




Fig. 5: Contracting two Seiherg-Witten half- geometries yields the 
Nekrasov- Okounkov partition function corresponding to a fermion flux p 
through the surface. 

The f/(l) Nckrasov-Okounkov partition function with fermion flux p (see flgure 5) is 
found by contracting the above fermion states 

Z^o(P;A) = (p|e^iex"-i|p). (5.31) 

The factors A can be pulled out of the exponentials by using the commutator [Lq, a±i\ — 
a±i. Up to an extra factor A~*'^/^ we find that 

Z^o(p;A) - (p|e^A2^°e^|p). (5.32) 

This has a nice geometrical explanation, since the left and right half geometries are con- 
nected by a tube of size A^ as in the factorized form of the complete U{1) geometry. The 
factor A^^° is the Hamiltonian that describes the propagation of the fermion field along 
the tube. There is no need to generahze this standard-CFT factor, since both patches are 
described by the same space-coordinate s. 

We also note that, as consistent with [4], the solution ^^(t; u) to PriP = equals the 
one-point- function 

(p-l|V'(t)|>V«) = J]r^'-"(p;it:„|>v«) ^r^'ei*-^ = (5-33) 

n 

where i?„ represents a Young tableau consisting of just one row of n boxes. 
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U{2) theory 

We apply now the above strategy for the U{2) geometry. We spht the corresponding 
curve into a left and a right half geometry, and for brevity focus just on the right part 
defined by 

Er: AH-^ ^ {v-U2){v-ui). (5.34) 

The 5-field canonically quantizes this equation into the second order differential equation 

PR^Pit) = {X\tdt-p2){tdt-p,) - AH-'}i;{s) = 0. (5.35) 

A change of variables z = 2t~^/^ followed by the ansatz = z~^P'^~^^^^ (l){z) and 

the rescaling z h- > {X/A)z transforms this differential equation into the familiar Bessel 
equation 

+ - I/' - z") 0(z) = 0, with i/2^(pi-p2)', (5.36) 

whose linearly independent solutions are given by modified Bessel functions and 
Ki,{z) of the first kind. The total solution in the original i-coordinate is therefore a linear 
combination of 



where p — p\ -\- P2- These modified Bessel functions have different asymptotics at infinity 
and relate to each other by going around the punctured disc C^. 
The second order differential operator Pr defines the "Dp^-module 

Mr^-^, (5.38) 

which we claim represents fermions on the quantum SU (2) Seiberg-Witten geometry. To 
check this statement, we have to find the fermionic state corresponding io Mr. So we 
asymptotically expand of the modified Bessel functions around t = in A: 

,/'2A\ ,i/4„^„/' 2A \ r jp-l) Ay/t , (^-l)(^-9) \h , 
k('^'^\ tV'^^^i f, , (ti-l) ^sft (,.-l)(/i-9) \H 1 



with ^ — Au 



2 
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Recall that equation (5.22) implies that any exponential function in the local coordi- 
nate — y/t near the puncture acts trivially on the vacuum state. Equivalently, this is 
true for any asymptotic series in y/i that assumes the value 1 at y/t — 0. In other words, 
we can forget about the complete expansion in ^/t\ Only the WKB pieces 

*'/^exp(±i^) (5.39) 

are relevant in writing down the fermionic state. This is exactly opposite to the matrix 
model examples, where the WKB-piece can be neglected and the perturbative series in A 
defines the fermionic state. 

The derivatives of the above sohitions have one term proportional to ip{s) (which we 
may forget about), and a term proportional to the derivative of the Bessel functions. The 
latter may be expanded as 

\X^/iJl 8 2A 2! -82 4A2 J 

aK(f) ^-^/^^rnP^^ ^^ I + (/.-l)(/x + l5) XH 

around Vi — 0. Again with the same reasoning only the WKB piece is necessary to write 
down the quantum state. Taking into account the extra factor in (5.37) the subspace 
is thus generated by the C(i)-module 

and blends (via the lexicographical ordening) into the fermionic state 

\y\;+) = v~P e^" (v^ A A A v"i A . . (5.41) 

on the cover. Here we used a cover coordinate v^^ obeying v~'^ = t, and rescaled the 
topological string coupling as A = A/2. is thus simply generated by a single function 

j/;^{v) = v'^e^" (5.42) 

Hence the fermions blend into the Bogoliubov state 

= (5.43) 

when p is an integer. 
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Note that the only modulus that appears in this expression is p. This represents 
the diagonal denoting the total fermion flux through the geometry. The moduli 

Pi and p2 measure the fermion flux through an internal cycle and are not visible in the 
result, because the final state sums over all internal momenta. In general any SU{2) 
Seiberg-Witten geometry with the same quantized p yields the same fermionic state. 

The fermionic (or dual) partition function is found by contracting the left and the 
right states, similarly as in the C/(l) example above. The left state is just the complex 
conjugate of the right one, so we find 

Z^o(p;A,A) = (p|eX"^eX"-^|p) - (ple^^A^^vK^ (5.44) 

The result is very similar to the U{1) example, up to the shift A i— > A/2. But notice 
that this fermionic state is written in terms of a single blended fermion. Decomposing 
this fermion into two twisted fermions makes it natural to insert an extra operator in the 
middle of the correlator, that measures the momenta of the two fermions through the 
^-cycles of the SW geometry. Weighting these momenta with a potential ^j, for i = 1, 2, 
yields 

Z^oi^„p;X,A) ~ (p|ei-e^«.A2^0ei"-b), (5.45) 

where H^^ = jJ2r ^(r+i/2) mod 2^r^lr = jiPi^i +^26)- This is the answer conjectured by 
Nekrasov and Okounkov in [55]. 

U{N) theory 

It is not difficult to extend this discussion to the U{N) theory (5.14), whose corre- 
sponding right half geometry we write as 

N 

Siv: A''t-' = l[{v-u,). (5.46) 

i=l 

Canonically quantizing this geometry and changing the coordinates z = {j}'^ t'^, brings 
us to the degree N differential equation 

It turns out that a solution to the above equation is given by a particular Meijer G- 
function, denoted G^^'^{z). The Meijer G-function is a complicated special function which 
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was introduced in order to unify a number of standard special function [58, 59, 60], and 
is defined in terms of a complex integral 



^ ' ) ~ 27ri A n,Wi r(i - 6, + 1) n,L„+, r(a, - 1) ^'-"'^ 



where L is a contour which goes from — ioo to +ioo and separates the poles of r(6j — i), 
for j = 1, . . . , m, from those of r(l — ai + t), for i = 1, . . . , n. 

It can be shown that the Meijer G-function solves the differential equation 

' q ^ \ 

W{zd, - bi) + {-ly-^-^+hHizd, - Uj + 1) G{z) = 0. (5.49) 
=1 j=i J 

So, indeed the Seiberg-Witten differential equation (5.47) is a special case of Meijer dif- 
ferential equation (5.49) with p — n — and q — N. Therefore the differential equation 
(5.47) is solved by 

i^i^)-GZ( ^ I^V (5-50) 

\ Pi,P2,...,PN ' y 

Similarly as before we claim that the P-module corresponding to U (N) Seiberg-Witten 
curve is generated by Pat. A subspace W corresponding to this P-module is thus generated 
by a solution iplt) and its derivatives in tdt- 

For p < q the Meyer differential equation (5.49) has a regular singularity at 2; = 
and an irregular one for z = 00. To extract the 1-branc fcrmionic state, we are interested 
in the behaviour around the irregular singularity, where t ^ 0. It turns out that one 
of the independent solutions of the Seiberg-Witten differential equation (5.47) has the 
asymptotic expansion [58, 59, 60] 

00 

A 1 — AT — ^ 

ip{v) e'^ v~ 2_^kjV~^ , (5.51) 

3=0 

around this singularity, which is conveniently written in the cover coordinate (— f)^ = 
= {j)'^ z. The other solutions are found by multiplying the coordinate v by A^-th 
roots of unity, and thus behave distinctly at infinity. As before, p = Ylf=iPi- 

To find the fermionic state corresponding to the U{N) Seiberg-Witten curve, we act 
with ip{v) on the Dirac vacuum. The positive power of v in the exponent of ijj{v) corre- 
sponds in the operator language to a-i, whereas lifts the Fermi level. The remaining 
series just contains negative powers of v which translate to a trivial action on the vacuum 
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in the operator formalism. Therefore, the above asymptotic solution and its derivatives 
(in tdt) blend into the state 

|>Vfl)=eKi|p), (552) 

with rescaled topological string couphng A = X/N. Like for the U{2) Seiberg-Witten 
geometry the dependence on the individual moduli Pi has dropped out. 

Similarly as in U{1) and U{2), in the present case we also find the U{N) Nekrasov- 
Okounkov dual partition function 

Z^ote; A,A) = (p|ei-e^e.A2^°ei"-». (5.53) 

This fermionic correlator is indeed the one postulated in [55]. For N = 1 or N = 2 the 
Meijer G-function specializes respectively to the exponent and Bessel functions, which 
reproduces the results derived in previous subsections. 

Although the normalizable moduli Pi disappear in the final I-brane partition function, 
they reappear when the state is unblended in terms of N single fermions 

as may be seen from (5.10) and (5.13). The charges Pi have an interpretation as the 
fermion fluxes through the tubes of the Seiberg-Witten geometry we started with. 

Actually, we find the same fermionic state when starting with any other Seiberg- 
Witten geometry whose fermion flux at inflnity is p. Hence one microstate in the total 
sum (5.54) can be interpreted as a fermion flux through an infinite set of geometries. This 
gives the state (5.54) as well as the partition function (5.8) the interpretation of a sum 
over geometries. 

5.3 Relation to topological string theory 

Nekrasov and Okounkov also derive a partition function for the 5-dimensional U (N) 
Seiberg-Witten theory compactified on the circle of circumference /3 [54, 55, 57] . It is 
given by a X-theoretic generalization of the 4-dimensional formula in equation (5.4). 

This 5-dimensional theory is closely related to the topological string theory by geo- 
metric engineering [61] on a toric Calabi-Yau background [62, 63]. Namely, the partition 
function of the topological string theory on an Ajv-singularity fibered over (whose toric 
diagram consists of A" — 1 meshes as in figure 6) is equal to the partition function of the 
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Fig. 6: On the left we see the five-dimensional U{2) Seiberg-Witten 
surface with fermion fluxes through its A-cycles, and on the right a cor- 
responding toric diagram. The fermion flux deforms the Kdhler lengths 
of the toric diagram as in equation (5.55). 

5-dimensional gauge theory given above, when the Kahler sizes of the internal legs are 

^ e'^(»'+-'^'), Qb-{^^Y^ (5-55) 

where Fi labels the vertical legs and B the horizontal ones. In the so-called gauge theory 
limit, when — 0, the topological string partition function reduces to the 4- dimensional 
Seiberg-Witten partition function. The corresponding B-model mirror geometry is of the 
form 

Xsw- xy-H{t,v)^0, (5.56) 

where H{t. v) = represents a Riemann surface of genus N — 1. In the gauge theory limit 
this surface becomes the Seiberg-Witten curve T^sw, parametrized as in the equation (5.3). 

In topological string theory it is natural as well to write down a dual partition function 
[4]. In a local B-model this allows the possibility of arbitrary fermion fluxes through the 
handles of the Riemann surface. In this setting it has been argued before that turning 
on a fermion flux is equivalent to deforming the geometry. More precisely, fermion flux 
parametrized hj V = PiBi changes the integral of the holomorphic 3- form over any linking 
3-cycle Ai, and thereby shifts the complex structure moduli Si — fl as 

Si + \pi (5.57) 

In the A-model fermion flux translates into wrapping D4 branes around 4-cycles, and 
thereby deforms the Kahler moduli. The I-brane partition function thus equals the dual 
topological string partition function. 
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Fig. 7: Three-cycles in the Seiberg-Witten U (2) -geometry. 

Because the Seiberg-Witten surface is embedded in C x C*, ^ and i3-cycles in the toric 
threefold will have topologies x S"^ and S^, respectively (they are drawn in figure 7). In 
particular, a basis of ^j-cycles can be chosen to reduce to the surface as the combination 
Ai — Ai^i. Now notice that the 3-cycle Ai with topology x S'^ is mirror to the vertical 
2-cycle Fi that connects the i-th and the i + 1-th horizontal leg. So turning on a fermion 
flux Pi through the z-th leg of the Seiberg-Witten geometry changes the complex structure 
parameter Si by an amount proportional to Oj — Oj+i. This explains the Kahler size Qp^ 
in (5.55) in terms of fermionic fluxes through the Seiberg-Witten curve, and in reverse 
why (5.54) may be interpreted as a sum over Seiberg-Witten geometries, or equivalently 
toric diagrams. So we conclude that the fermionic interpretation in 4d of Nekrasov and 
Okounkov is dual in 6d to the fermionic interpretation of the topological string, and has 
a deeper interpretation in terms of P-modules. 

Five- dimensional U{1) theory 

Quantizing a five- dimensional Seiberg-Witten geometry yields a difference (instead of 
differential) equation. Working out I^-modules for these geometries we leave for future 
work. Let us treat one example in detail though. The five-dimensional right U{1) Seiberg- 
Witten half-geometry 

S^'^ : pAe-f^H-' + e-^" -1 = (5.58) 

may be drawn as a pair of pants. In the field theory limit /5 — it reduces to the familiar 
equation At~^ = v for the right-half Seiberg-Witten geometry (with m = 0). 

In the B-model the most general state assigned to a local pair of pants geometry is 
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given by a Bogoliubov state [4] 



oo 



l>V>=exp[j] Yl «^nV'i„.-l/2V'jl_V2]|0>, 



(5.59) 



i,j m,n=0 



where the index i — 1,2,3 describes the fermion field on the three asymptotic regions of 
the pair of pants, and the coefficients are determined by a comparison with the A-model 
topological vertex. This exponent can be expanded as a sum over states (see figure 8) 



where the fermion fiux is conserved: pi + P2 + Ps — 0. To describe the 5d Seiberg-Witten 
U{1) geometry we won't need this state in full generality. 

Ri 



Fig. 8: The B-model vertex (on the left) may be expanded as a sum 
over fermionic states \pi,Ri) \p2, R2) ® |P3, Rz), with pi +P2 +Pz = 0; 

corresponding to a conserved fermion flux through the pair of pants. The 
five- dimensional right-half Seiberg-Witten geometry (on the right) with 
charge p only has one partition Rj^Q. 



The B-field quantizes this geometry into the difference equation 




(5.60) 




R 



P 




P(t)^(t) = (/^Ae-'^V^ + e' 



,l3Mdt 



1) ^{t) = 0. 



(5.61) 



Its fundamental solution is the quantum dilogarithm 




(5.62) 



As an intermezzo, notice that quantizing the equation 



log (1 - pAe-f^H-^) , 



(5.63) 
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which is just a rewriting of equation (5.58) for E|^, we find a differential equation which 
may be interpreted as the WKB approximation of difference equation (5.61). A funda- 
mental solution of the differential equation is given by the genus disc amplitude 

n>0 

Acting with the five-dimensional dilogarithm on the Dirac vacuum state yields the 
fermionic state 

l>^)"a) = -pEifz5STl0)- (5-65) 

n>0 ^ ^ 

This describes a subset of |yV) where only the quantum number Ri is non-trivial. Sum- 
ming over all external states of the form 

I ® ® |0,«), (5.66) 

incorporates a fermion flux p through the pair of pants. In the field theory limit /? ^ 
the resulting state reduces to the familiar four-dimensional state 

exp(Q;_i/A)|p) ® \p, •) ® |0, •). 

The partition function is found as the contraction of the left and right 5d half- 
geometries. (Or equivalently in the topological B-model by inserting a propagator [4].) 
This yields the fermionic correlator 

(0|f+f_|0) = (0|r+(/?A)2^°r_|0), (5.67) 

with 

~ ^ (/9A)l"la„ s-^ an 

r^-^-P E H(i-e^An) r,-exp^ (5.68) 

±n>0 ' ^ ±n>0 ^ 

Indeed, the result equals the five-dimensional U{1) partition function 
that was found by Nekrasov and Okounkov in 
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6 Discussion 



In this paper we argued that the fundamental objects underlying various systems in 
theoretical physics are chiral fermions living on quantum curves. In our formulation the 
quantum curve is defined, similarly to an affine classical curve, in terms of an equation of 
the form P{z, w) = 0. Its crucial feature, however, is the non-commutative character of 
the coordinates z, w. These quantum (or non-commutative) curves generalize the classical 
curves that come up in the standard formulation of a given topic. Examples of such 
classical curves can be found in the theory of random matrices, c = 1 string theory, 
Seiberg-Witten theory, and more generally in topological string theory. Semi-classically 
their (genus one) free energy is computed as a fcrmionic determinant on the classical curve. 
In our approach chiral fermions on the quantum curve generate the all-genus expansion 
of the free energy with respect to the non-commutativity parameter A. 

Fermions on a non-commutative curve can be realized physically within string theory 
as masslcss states of open strings on an intersecting brane configuration in the presence 
of a B-ficld. This idea was already put forward in [5]. In this paper we have exploited 
this I-brane system in a few important examples. 

First of all we showed, reinterpreting the results in [38], that I-branes and "D-modules 
provide an insightful formulation of matrix models. This quite general statement is also 
appealing when certain matrix model limits are considered, such as a double scaling limits. 
In this case one recovers an I-brane formulation of minimal string theory and topological 
gravity. Secondly, we discussed how to reformulate c = 1 string theory in the framework 
of D-modules. 

Finally, we discussed supersymmetric gauge theories. Using P-module formalism we 
derived fermionic expressions for the partition function of the M = 2 gauge theory, re- 
producing the dual all-genus partition function introduced in [55]. We considered mainly 
4-dimensional Seiberg-Witten geometries with unitary gauge groups, and explained only 
the simplest U{1) example in the 5-dimensional theory. It would be insightful to extend 
these results to other gauge groups and include matter content. It is clear that this 
should be possible, as these aspects of the 5-dimensional Seiberg-Witten theory are 
captured by topological string theory on toric manifolds. The latter system can be solved 
in fermionic B-model formulation of the topological vertex [4] which is equivalent to the 
I-brane fermions [5]. Nonetheless, finding the quantum I-brane curve representing such 
configurations appears to be a nontrivial task. 
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In all these examples we were able write down a "D- module that, through the pre- 
scription in section 2, yields the all-genus partition function. Especially the matrix model 
examples made it clear that this "D-module can be quite non-trivial in general. Only for 
the simplest curves, such as those appearing in double scaled matrix models, the "D-module 
can be found by canonically quantizing the classical spectral curve. 

In the process of unraveling the P-module structure in both sets of examples, we 
noticed some crucial differences. While the WKB piece of the P-module generator can 
be ignored in finding the all-genus matrix model partition function, we discovered that it 
plays an eminent role for the Seiberg-Witten geometries. Another distinction is the differ- 
ence in (non-)normalizable modes. While the potential W parametrizes non-normalizable 
modes that appear in the X>-module as parameters, in contrast, the normalizable modes 
in the Seiberg-Witten geometries are eaten by the D-module, and only visible as a sum 
over internal fermion fluxes in the geometry. On the other hand, varying the "D-module 
with respect to the non-normalizable modes yields differential equations which relate to 
isomonodromy and the Stokes phenomenon. 

Even with this rather broad set of examples, a few major questions remain. First 
of all, we cannot give a recipe in general how to find the quantum curve underlying 
a certain problem. Secondly, it is not obvious that our prescription is independent of 
the chosen parametrization of the classical curve. As we noted in the example of the 
classical curve zw — 1, different parametrizations can lead to different quantum curves 
that nonetheless yield the same partition function. This should hold in general cases as 
well, as topological string theory associates a unique all-genus partition function to a 
given curve. Thirdly, we haven't exploited some of the advantages of using P-modules 
instead of differential equations. One of the main advantages is some independence on 
the way the differential equation is written down. It would be very interesting to try to 
match this freedom with the choice of parametrization for the classical curve. Finally, we 
have only discussed examples with one or two local patches. It would be highly insightful 
to study more general examples. 

While in this paper our focus has been to associate a A-perturbative quantum state 
to a spectral curve, we noticed that "D-modules in fact contain non-perturbativc informa- 
tion. These bits get lost when we turn the D-module in a fermionic state by making an 
asymptotic expansion of the P-module generators in A. This is in line with the discussion 
in [27] , where it is argued that non-perturbative effects drastically modify the non-trivial 
target space curve into a complex plane. Non-perturbative effects in matrix models, as 



55 



well as in the topological string theory, were also recently discussed in [64, 65]. Especially 
interesting in this respect is [66] , where a non-perturbative partition function is proposed 
that is very similar to the I-brane partition function (2.12). 

In the step where we turn a X>-module in a quantum state, a choice of boundary 
conditions has to be made. This implies that final states are troubled by the Stokes 
effect: solutions that decay faster can be added at no cost and the state changes when 
one crosses certain lines in the moduli space. This suggests that the P-modules we 
studied in this paper may help in the understanding of wall-crossing phenomena in the 
corresponding M — 2 theories [67, 68] . 

More mathematically, our formalism is deeply connected with quantum integrable 
systems and the geometric Langlands program [69, 8, 70, 21, 22, 71, 72], while approaching 
these topics from a string theoretic perspective. Especially interesting in this respect is 
our quantitative approach, that allows us to associate quantum invariants to spectral 
curves. In the future we hope to make this link even more concrete. 

Specifically, it would be enlightening to have a better description of the non- 
commutative fermionic CFT on a given quantum curve. It is interesting to find out 
whether this relates to the WZW models based on opers in the geometric Langlands 
program: as so-called Hecke eigensheafs these generate examples of the Langlands cor- 
respondence. And, to discover the relation with the interacting bosonic CFT's that give 
another perspective on these intersecting brane configurations [73, 74, 75] as well as [76]. 
In particular, both models determine a set of recursion relations. It would be helpful to 
compare them. 

A clear physical realization of quantum curves and the associated well-defined math- 
ematical formulation in terms of D-modules are great advantages of our approach. In 
consequence it can be applied to mimerous situations mentioned above and yields definite 
quantitative results. Nonetheless, the idea of quantum curves is not new and earlier at- 
tempts of their formulation appeared before in physics and mathematics. It is worthwhile 
to recall how those attempts relate to our formalism. 

The notion of quantum or non-commutative geometry has also been introduced by A. 
Connes [77] . His approach relies on replacing the algebra of functions on a manifold by 
a non-commutative C*-algebra. In this context a program of developing a theory of non- 
commutative Riemann surfaces, from the point of view of geometric quantization [78], 
was advanced in [79]. Independently of this program, also some particular examples of 
low genus non-commutative Riemann surfaces have been analyzed in literature. In genus 



56 



zero they include the so-called Podles sphere [80] and more generally fuzzy spheres [81], 
which also found vast application in string theory. In genus one, one can consider a non- 
commutative torus which arises naturally in a certain realization of M-theory known as 
Matrix theory [82, 83] . Non-trivial S-field is an essential ingredient in a realization of these 
systems. It would be interesting to see if they could be related to I-brane configurations. 
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A Infinite dimensional Grassmannian 

In this section we introduce an infinite dimensional Grassmannian and its description 
in terms of the second quantized fermion field (we learned this material e.g. from [3, 13, 
14, 15, 16]). 

Grassmannian and second quantized fermions 

The space H — C{{z~^)) of all formal Laurent series in can be given an interpre- 
tation of a Hilbert space. Basis vectors 2;", for n e Z, correspond to one particle states of 
energy n associated to the Hamiltonian zd^. This Hilbert space has a decomposition 

H = n+®n-, (A.i) 

such that the first factor H+ — C[z] is a subspace generated by 2;°, z^, z^, . . ., while H- 
is generated by negative powers z~^, z""^, Consider now a subspace W oi H with a 
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basis {wk{z)}ken- We say it is comparable to Ti.+ , if in the projection onto positive and 
negative modes 

= ^{w+)ijz^ + ^{w-)ijZ-^ (A.2) 

i>o j>o 

the matrix w+ is invertiblc. The Grassmannian Grg is the set of all subspaces W C C((^)) 
which arc comparable to 7i^. 

In what follows wc take much advantage of the correspondence between Gtq and the 
charge zero sector of the second quantized fermion Fock space J-'q. In this correspondence 
the subspace 7i+ is quantized as the Dirac vacuum 

|0) A^^ A^^ A..., (A.3) 

with all positive energy states filled. The fcrmionic state associated to the subspace W 
with basis wo{z), wi{z), W2{z), ... is represented by the semi-infinite wedge^ 

\yV) ^ Wo Awi AW2 A . . . (A.4) 

which is an element of the fiber of a determinant line bundle over the element W G Gr 
(and therefore determined up a complex scalar c). 

To make contact with the usual formulation of the second quantized fermion Fock 
space, we can identify the differentiation and wedging operators with the fermionic modes 

^n+l = ^ C+i = A . (A.5) 

These half-integer modes are annihilation and creation operators which arise from a de- 
composition of the fermion field ip{z) and its conjugate ip*{z) 

i;{z) = Yl ^*(^) = H rrZ-^''^, (A.6) 

rez+i 

and they obey the anti-commutation relations {■0r,'0-s} — ^r,s- 

For subspaces W G Gro the determinant of the projection onto H+ is well defined and 
can be expressed as 

detw+ = (0|>V). (A.7) 

More generally, one can consider the Fock space which sphts into subspaces of 
charge p 

^ = ^p. (A.8) 
^Actually, we have to tensor with to make the state fermionic. 
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Each subspace J^^ is built by acting with creation and annihilation operators on a vacuum 

|p) = z^ K zP+^ A zP+^ A..., (A.9) 



with the property 

ipr\p) — for r > p, 

ip;\p) = forr>-p. (A. 10) 

The Fermi level of the vacuum \p) is shifted by p units with respect to the Dirac vacuum 
|0). This fermion charge is measured by the U{1) current 

J{z) =: i/j{z)r{z) := 5]a„^-"-\ (A.ll) 

n 

whose components = : ipr'4'n-r satisfy the bosonic commutation relations 

[am, a-n] = rn6m,n- (A.12) 

With each subspace W C C((^)) comparable to the one generated by {z'^)k>p one can 
associate a state | W) G of charge p. This charge is equal to the index of the projection 
operator pr+ : W H+. 




P 



Fig. 9: Elements of the Fock space T are in a bijective correspondence 
with Maya diagrams. The bottom line represent a Maya diagram corre- 
sponding to a fermionic state with charge p. As illustrated it is char- 
acterized by a two-dimensional partitions R located at position p. We 
therefore denote the state as \p, R) G J^. 

A state in the Fock space has also a simple representation in terms of the so- 
called Maya diagram (see figure 9). Black boxes in such a diagram represent excitations, 
whereas white boxes are gaps in the energy spectrum of the fermion. The charge of a 
state is given by the number of excitations minus the number of gaps. Fermionic states or 
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Maya diagrams of a fixed charge p can also be associated to two-dimensional partitions. 
In particular in p = sector the state 

1=1 

corresponds to the partition R — . . . ,Ri) such that 

ai = Ri-i, hi = R\-i. (A.14) 

In what follows a state corresponding to a partition R of charge p is denoted as |p, R) . 

Flow on the Grassmannian 

There is an action on the Grassmannian defined by multiplying a basis vector Wk{z) 
of W by a power series f{z) — J2 fn^^ that vanishes at 2; = 0. 

f{z)\W) =^WoA...(^fc-i A/-WfcAWfe+i.... (A.15) 

k 

When we write Wk{z) in terms of the basis {z^)i^z this action is encoded by the multipli- 
cation by an infinite matrix in glao-, whose {i,jy^ entry is given by On the fcrmionic 
state |VV) a multiplication by 2;" translates into a commutator with the bosonic mode q;„, 
since increases the fermionic mode number by 

[an,i>r] = i>r+n- (A. 16) 

Multiplication by a power series f{z) therefore translates to the operator 

/ = Zl-^'^t"'^'*] ^9loo (A. 17) 

n 

on the Fock space. 

Exponentiating the action of gl^o yields the group G/qo- An element g{z) = exp{f{z)) 
of this group acts on \yV) by multiplying all its basis vectors 

g{z)\W)^g-WoA...Ag-WkA.... (A.18) 

From the fermionic point of view this action is given by conjugating each basis vector Wk 
with the element 

g = exp (^^/nOn) = exp (^^ dz f{z)J{z)^ e Gl^o- (A.19) 



60 



We call r the group of exponentials g{z) : ^ C*. An important subgroup of F is 
the group r+ of functions go : ^ C* that extend over the disk Dq — {z : \z\ < 1}: 

r+ = {go -.Do^C*: go{0) = 1}. (A.20) 

Another subgroup is the group r_ of functions g^o ■ ^ C* that extend over the disk 
-Doo = {-2 e C U {oo} : 1^1 < 1}: 

r_ = {g^ : Doo ^ C* : ^00(00) - 1}. (A.21) 

Any (7 e r can be written as an exponential exp(/). When g eT^ the function / vanishes 
at 2; = 0, and when gi e r_ it vanishes at 2; = 00. 

r+ and r_ have different properties when acting on Grassmannian. The action of 
r_ is free, since any W G Gr has only a finite number of excitations. On the contrary, 
r+ acts trivially on a vacuum state \p). Although the action of the groups r_|_ and r_ 
on a subspace W is commutative, as it is just given by multiplication, as operators on 
the fermionic state |>V) it matters which element is apphed first. This introduces normal 
ordering ambiguities. 

An element 

g{t, z) = exp I J] hz'' 1 = exp {f{t, z)) e r+, (A.22) 
\fc>i / 

defines a linear flow over the Grassmannian Gr. On the Fock space it acts as an evolution 
operator 

U{t) = exp ^J{t, z)J{z)^ . (A.23) 

The determinant det(yV)+ is not equivariant with respect to the action of r_|_. The 
difference is measured by the so-called tau-function 

det jg-'wu m{t)m 

^^^^^ = ^-det^^ = g-HOm ' ^^-'^^ 
which yields a holomorphic function r : r+ — > C. This can be regarded as a wave function 

of \yv). 

Blending 

So far we considered the Hilbert space TC = Ti.'^^^ of functions with values in C. More 
generally, one can consider a Hilbert space H^"'^ of functions with values in C". Let 
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(Q)i=i,...,n denote a basis of C". For each n there is an isomorphism between 7Y^"^ and H 
given by the lexicographical identification of the basis 

€iZ^ ^ (A.25) 

This isomorphism is called blending. 

In the fermionic language the Hilbert space hfts to the Fock space of n fermions 
i — 1, . . . ,n, each one with the expansion (A. 6) and such that 

{i^l^\rs^^^S,jSr,-s. (A.26) 

Now blending translates to the following redefinitions of these n fermions into a single 
fermion ijj 

where 

2i — n — 1 / » 

ft = (A.28) 

Blending can also be expressed in terms of two-dimensional partitions introduced 
above. Consider n partitions R(i) of charges pi, with YliiPi = Pi corresponding to states 
in n independent Hilbert spaces of fermions ■0'^'^. Associating with each such partition a 
state of a chiral fermion we have a decomposition 

n 

|p,R) = (g)b„i?(,)), (A.29) 

i=l 

and the blended partition R of charge p, corresponding to a state in the Hilbert space of 
the blended fermion ^, is defined as 

{n{pi + - m) + i - 1 I m e N} = {p + Ri^ - X I X e N}. (A.30) 



B Some background on P-modules 

The theory of "D-modules was introduced and developed, among others, by I. Bernstein, 
M. Kashiwara, T. Kawai and M. Sato, to study hnear partial differential equations from an 
algebraic perspective [9, 10, 11, 12]. Currently this is a very active field, with connections 
and applications to many other branches of mathematics. 
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"D-modules are defined as modules for the algebra of differential operators V. In 
general, in a local C" patch with complex coordinates {zi, . . . , Zn), the operators Zi and 
^z^ represent the n*^ Weyl algebra. The operators P e T> are of the form 

il,...,in 

With a set of operators Pi, ... , Pm € one can associate a system of differential equations 

Pi* = . . . = Pm^ = 0, (B.2) 

where ^ takes values in some function space V. An algebraic description of solutions to 
such a system can be given in terms of a X>-module Ai determined by an ideal generated 
by Pi,...,P^eP 

M = r. (B.3) 

V-{P,,...,P^) ^ ' 

The advantage of considering such a D-module is, firstly, that it captures the solutions to 
the above system of differential equations independently of the form in which this system 
is written. Secondly, it is also independent of the function space V - be it the space of 
square-integrable functions, the space of distributions, the space of holomorphic functions, 
etc. 

Nonetheless, having chosen a particular space V one is interested in, the space of 
solutions is simply given by the algebra homomorphism 

Hom^,(M,V). (B.4) 

E.g. holomorphic solutions to the differential equation P^{z) — can be captured as a 
homomorphism of "D-modules 

= ^ - (^c, (B.5) 

with Oc the algebra of holomorphic functions on the complex plane C. Indeed, define a 
map that sends the element 

[1] e M ^ e Oc- (B.6) 

This is well-defined because every element P' G DP is mapped to zero (remember that ^ 
fulfills P* = 0), and it is a bijection; conversely, any map M. to Oc is determined by a 
holomorphic solution to the differential equation P^! = 0. 
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An important notion is a dimension of a "D-module. The so-called Bernstein inequality 
asserts that a non-zero "D-module M. over the n*'* Weyl algebra has a dimension 2n > 
dimA^ >n. In particular, V considered itself as a "D-module has a dimension 2n. On the 
other hand, dimC[xi, . . . , Xn] — n. For a non-zero P eT>, dimD/T>P — 2n — 1. 

A special role in the theory of "D-modules is played by the so-called holonomic T>- 
modules, which by definition have a minimal dimension n. In particular they are cyclic, 
which means of the form {D'i/ : D e V}, i.e. they are determined by a single element 
^ E A4 called a generator. 

In the context of the I-brane in we are just interested in the 1** Weyl algebra 
(2.15) of dimension 2. In this case we immediately conclude that the module V/VP has a 
dimension n = 1 for any non-zero P, and is thus holonomic and cyclic. It can be realized 
as 

M = -.DeV}, (B.7) 
where the generator ^ is a solution to the differential equation P^ — 0. 

Flat connections 

More generally, P-modules are defined as differential sheaves on any variety X. The 
sections of the sheaf T>x over an open neighbourhood U are given by linear differential 
operators on U. Therefore, both the structure sheaf Ox (of holomorphic functions) as 
well as the tangent sheaf Tx (whose local sections are vector fields) may be embedded in 
Vx 

Ox-^Vx^ Tx. (B.8) 

In fact, Vx is generated by these inclusions. 

A sheaf on X is defined to be a left module for Vx when i> • s e A4, for any v G Vx 
and s e A4. Furthermore, it has to fulfill 

v-{fs)^v{f)s + f{vs) (B.9) 
[v ^iW] ■ S — V ■ {w • s) — w • {v ■ s) 

for any v e Vx, f e Ox and s E A4. Suppose that At is a left "Dx-module whose sections 
are the local sections of some vector bundle V (this encomprises all Px-modules that are 
finitely generated as Ox-modules) . Then the action of Vx defines a connection on V as 

V^{s) = vs, (B.IO) 
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whose curvature is zero. So a "D-module structure on the sheaf of sections of a vector 
bundle V defines a fiat connection on this vector bundle. And conversely, any module 
consisting of sections of a vector bundle V with fiat connection Va, has an interpretation 
as a "D-module defined through the action of the flat connection. Therefore, a "D-module 
is in general just a system of linear differential equations, changing from patch to patch 
on X. This is known as a local system. In the main part of this paper X is just C or C*. 

C Relation to quantum integrable systems 

In this article we focus on smooth curves that are given by an equation of the form 

E: H{z,w)^w'' + Un-i{z)w''-^ + ... + Uo{z)^0, (C.l) 

where 2; e C (or C*) and w & C These play a prominent role in integrable systems as 
spectral curves. It is a degree n cover over C (or C*) 

E C T*C 

i^r (C.2) 

C 

with possible branch points (from now on we restrict to z E C for simplicity in notation). 
The spectral curve is embedded in and equipped with the (meromorphic) 1-form 

T) = \wdz\j:. (C.3) 
A 

Our notion of a quantum curve agrees with a notion of quantum spectral curves in 
this context. Let us say a few words about this. 

1/2 

Fermions on E transform as holomorphic sections of a line bundle L (g) Kj! , provided 
by the D6-brane. The pair (L, 77) on E pushes forward to a couple 

TT* : {L,r]) {V = ■K^L,(p = TT^rj) (C.4) 

on C under the projection map tt : E — > C. So y is a rank n vector bundle on C, whereas 
is a holomorphic 1-form valued in gl{n).^ Such an object is called a Higgs field. It 
endows V with the structure of a Higgs bundle. Setting the characteristic polynomial 

det(r; - (t){z)) = (C.5) 
^In other words, 4> G H°{C,EndV O Kc). 
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returns the equation for the spectral curve. The push-forward map tt* sets up a bijection 
between spectral data and (stable) Higgs pairs 

(E,L) ^ iVA). (C.6) 

The moduh space of stable Higgs pairs is an algebraically completely integrable system, 
known as the Hitchin integrable system 

A 'D;^-module (as in [20]) corresponds to a A-connection 

V^ = Xd,-A{z), (C.7) 

which is defined through the Leibnitz rule V\{fs) — fV\{s) -\-\s®df for any function / 
and section s. 

Semi-classically, such a A-connection Va reduces to a 1-form Vo(-z) with values in gl{n) 

Va^Vo, (A^O). (C.8) 

We just encountered this object as a Higgs field 0. Moreover, we explained with (C.4) 
that a Higgs (V, (f)) and spectral data (E, L) provide equivalent information. In particular, 
the spectral curve can be recovered by the determinant of the Higgs field. This implies 
that A-connections quantize spectral data.^ 

It tells us exactly which requirements a "D-module quantizing the I-brane configuration 
has to satisfy. Fermions on a degree n spectral curve have to transform under a rank n 
A-connection Va on C, whose semi-classical A ^ limit is given by the Higgs field 

Vo = 7r,(7?). (C.9) 

A simple example of a A-connection is given by 

V ^Xd,- A{z), (C.IO) 

with A{z) = 71^,(1]). Its determinant is a degree n differential equation that canonically 
quantizes the defining equation for S. 
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